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ABSTRACT 


S  me  Nonlinear  Feedback  Control  Systems  under  high  gain  Gon- 
dii,  s  exhibit  the  phenomen  n  of  subharmonic  instability s  or  contain 
subharn  ,i.s  .  f  the  fundamental   Jtput  frequency,  A  general  discussion 
of  subharmonic s  in  nonlinear  systems  is  followed  by  an  investigation 
of  ripple  instability  Ln  a  forced  limit -cycling  voltage  regulator  con- 
taining  a  thyristi  r  <  —  SCE  briJge  utilizing  an  ON -OFF"  switching 
scr eme . 

A  digital  computer  pr  ^r^m  is  used  to  simulate  the  dynamic 
response  of  the  system  under  different  loading  conditions  and  for 
liffere  t  reference  voltage  levels,, 
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1.      Introducti 

There   are  many  physical    applications  where   a  very   stable  ditect 
current   power   supply   is  desirable.      This    is    usually  accomplished  by 
I  -king    -x  multiphase    .  lr    ■:      1 1  ig   curre/,    signal,    full  wave    rectifying 
it   and   then  passing    it   through  what    is  essentially  a   low  pass    filter. 
The   output   then   is    some  DC    level   on  which   rides   a   small  value    of 
"ripple".      The    ripple    frequency   is   some  multiple   of   the    source    input 
signal    dod    is   due    to  the    filter's  being  non-ideal. 

For  m   st    applications   the   amplitude   of   the    ripple   voltage    is 
not   so  great   that    the   power   supply   cannot  be   used  and  considered 
t      be   a   stable  DC    source, 

H   wever,  with   the   advent   of    the   c ■  impute r  and   switches   that 
need  only  millivolts   of   fluctuations    to   switch,    it    is   desirable    to 
have    power    supplies   tb^.t   are   more    tightly  regulated  and  effect   a 
East    recovery   fr  >m   fluctuations    in  source    supply  and  noise   pertur- 
bations „     This    scheme    is    usually  accomplished  with  higher   loop  gain. 

When  the    loop  gain   increases    to  higher  values,   however,    it 
has  been  witnessed   that   the   amplitude    of   the   ripple    increases   and 
c    ntains   components    of    frequencies   that   are    lower   than   those  wit- 
nessed under   lower  gain  conditions. 

F    wer   supplies   that   employ   silicon  controlled   rectifiers 

(thyristors)   which   are   of   the   phase -control   type   are   known  to  contain 

J  L 
these    subharmoriics   under  high   gain  conditions. 

The   nature    of    this   piper    is    to    investigate   and  verify   the 
occurrence    _f   subharai  nics    in   a   f   reed   limit-cycling    regulator  em- 
ptying   silicon   c    ntr<  Lied   rectifiers    that    are   not    phase    fired. 
Wliile  verifying   the    fact    th«.t    they  do  exist  ,    an    investigation   is 


m  ide  in  an  attempt  to  determine  how  these  subharmonics  are  caused , 
what  frequencies  are  presents  and  under  what  conditions  they  occur„ 

The  ultimate  gi  il  •  f  such  Investigation,,  of  course ,  would  be 
C>  di  sever  means  to  pre  i  k  r.  -J  -:      t<    .r  ranee  of  subharmonics  }    their 
frequencies  and  amplitudes,,  with  hopes  of  designing  means  to  com- 
pensate the  system  so  that  these  subharmonics  do  not  present  too 
great   f  a  degrading  effect  >n  the  output  of  the  regulated  power 
supply o 

2„  A  general  study  of  subharmonics  in  a  nonlinear  system. 

If  a  linear  system  is  forced  by  a  sinusoidal  input,  the 
resulting  steady -state  output  is  also  a  sinusoid.   However,  a  non- 
linear system  driven  by  a  sinufpid  does  not  produce  a  sinusoidal 
<  utputs  but  the   jc^j!  c  nt-niis  harmonics  which  are  frequencies  that 
are  integral  multiples   f  the  driving  frequency. 

In  some  cases  the  output  may  contain  one  or  more  frequencies 
that  ire  I  wer  than  the  forcing  frequency.  These  frequencies  are 
characteristically  integral  submultiples  «  f  the  driving  function 
and  -.re  called  subharmonics.  The  subharm  .  ics  may  be  small  in 
amplitude  when  compared  to  the  driving  frequency  or  may  be  so  great 
i  i  c  Mirparison  that  the  forcing  frequency  may  be  neglected. 

Unf  >rtunatelys  there  are  no  general  rules  available  by  which 
one  may  ascertain  whether  the  occurrence  of  subharmonic  oscillations 
is  possible.   Many  references  on  the  subject  suggest  that  the  question 
of  predicting  subharmonics  can  be  answered  by  c  insidering  whether 
or  not    o  lus  Ldally  f  rued  differential  equation  of  the  form 


J-  +  /<  Cr   y  -f(Lr)  -  3  Cos  pt 


(2-1) 
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can  I  ive    ^    &.  Luti<        it     i    frequency     >f     l)  /N ,   where  N    is   an   integer. 

Equ   'i     i    (2=1)    is  a  classical  equation  of  the   theory  of  non- 
linear  systems  and    is   kn  wn   as  Duffing' s  Equation. 

It  h.4us  beer,    found    from  experimentation  that   the   occurrence   of 
subharmonic    oscillations   in  physical    systems    is   strongly  dependent 
on  the    starting   conditions.      The   amplitude   and   frequency  of  the 
driving    force  must    fall  within  certain   definite    limits.      Because   of 
this    st  r    :ig  dependency  on    initial    conditions    it  would   appear   that 
sub-harmonics  would   only  be    a    trs    sient   problem.     However,    in   a 
n     ilinear   system  the   oscillation    is   non- sinusoidal  and  contains 
harmonics   of   the    fundamental    frequency.      It    is    possible   to  maintain 
the   oscillation   in  a   steady-state 9   under   s^me   conditions ,   by 
supplying  energy    t      the   system  at   any  of   these  harmonic    frequencies 0 
Hences    the   driving    force    is   a  harmonic    of   the    fundamental    frequency;, 
or  the    fundamental   frequency   is   an   integral   submultiple   of   the 
forcing    frequency.      This    then    is   a   condition    for   subharmonic 
gene rat  ion. 

2.1     Buffing's  equation 

Perhaps   the  most  exhaustive    study  of   the    subharmonic   phenomenon 
has  been   accredited  tn  Hayashi.        The    following  examples   of   subhar- 
m  nic     »scillation  applying  Duffing5 s  Equati     «  have  been   taken   from 
h  i  s  w<  r  k  „ 

Goi  sidei    the    fundamental  equati. 

ff    +  **    jg  *  A"«->  ="5CoS  ri  t,     (J.l.i.4'    )        (2.2) 

in  which  2  <$   is  a  c  astant  damping  coefficient  and  f(v)  is  a  term 
representing  the  nonlinear  r^st  ring  force. 
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The  period  of  the  forcing  function  is  2  ft'  /i)5  and  the  subhar- 
monic  L/|>  has  a  period  2 1Y  and  may  be  expressed  by  a  linear  com- 
bination of  sin  V   and  cos  t"  „ 

What  follows  is  an  investigation  of  the  relationship  between 
the  nonlinear  characteristic  expressed  by  the  term  f(v)  and  the 
orde r  1/0  of  the  s ib h -> r m<  nic  oscillati on s . 

Consider  the  restoring  force  to  be  the  polynomial 

^)r£  C^m  c»<r+*S+'s<r\  ■■■  (2.3) 

where  c.s  c.5  c„soo.  are  constants  which  are  determined  by  the  non- 
linear characteristics,,   These  constants  are  subject  to  the  constraint 
that 

c.   +  c,     +  c  -  +  e,   +  „  „ ,  =   1  (2-4) 

In  steady   state    one  may  assume   that    the    solution  of  equation 
(2-2)   has   the    farm 

y^     r    £    +    /.'sen    £   +   ucost-h    urease  (T 
Ihe   approximati  <n    f   r    uT    is  given  by 

U)"    =  1-    pi-         °  (2-6) 

which  is  legitimate  as  long  as  the  nonlinearity  is  small. 

By  substituting  (2-5)  Lnt   (2-2)  and  equating  the  coefficients 
of  the  sintl   and  cost*   terms  separately  to  zero,  the  following 
results  are  obtained  ac  :  rding  tu  the  form  of  the  nonlinear  character- 
i  s  t  ic s  f  ound  i * ,  ( 2 - 3) . 


12 


3 
Case  1.    f(v)  =  c_v  +  c  ,v  „ 

The  non linearity  is  symmetricals  or  f(v)  is  an  odd  function, 
and  in  this  case  the  constant  term  z  can  usually  be  disregarded „ 
Then  f<  t     =  2s4s5.o.:  §ubscir  ition  as  menti  ne  J  above  yields 
C  X  -  3/4  C  7*4  ^  *)  -  3/a  UT^*  +   k  U  £  o 

[i  -  V^  (  y"N  >,'j  ~  V*  "r*J*   -  ktj««  (2.7) 

where   K  r  '  riCa   .   By  simultaneous  solution  one  obtains  that 

k(x2  +  y2)  =  0  (2-8) 

Equation  (2-8)  implies  that  the  amplitude  of  subharmonic 
oscillations  is  zero  as  long  as  k#03  or  damping  is  present.  Hence 
subh-rnK  oiic  s  of  order  I./2,  1/4 s  1/5  s .  „  „  cannot  occur  in  this  case. 
But  real  roots  of  x  and  y  which  do  not  simultaneously  gc  to  zero 

may  be  obtained  for  0    -   3.   Thereby,,  one  may  conclude  that  sub- 

3 
harmonics  of  order  1/3  may  occur  when  the  nonlinear  term  c _ v   is 

present. 

Case  2.    f(v)  =  c.v  +  c,.v   +  c  v  . 

Since  the  n  all  .-.  rity  is  unsymmetrical  the  constant  term  z 
in  (2-5)  must  be  c  nsidered.  It  is  found  under  rigorous  solution 
that  subharmonics  of  the  order  1/2  may  exist  in  this  case. 

Case  3.    f(v)  =  c_v  +  :  _  v  . 

3 
Airi-     ugh   the   cubic    term  c  ,v     d  <es  n<  t    appear    in  the  nonlinearity, 

detailed    Lnvestigati    n   reveals   that   subharmonics   of   order    1/3   still 

a  ,r     a]     ig  with   those      f   order    1/5. 
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After  establishing  theoretically  the  possible  occurrence  of 

subharmonic  solutions  to  Duffing" s  Equation,  Hayashi  presents  the 
results  of  several  experimental  investigations ,  a  synopsis  of  which 

f  ILL  ws . 

2.1.1  Experimental  investigation  of  subharmonic  of  order  1/3. 

The  experimental  data  was  obtained  using  an  electric  oscil- 
latory circuit  c  citaining  a  saturable -core  inductor  and  a  capacitor 
in  series.  The  circuit  equ- r 1  a    takes  the  form  of  Equation  (2-2) 
■'  an  a  60  cps  a.c.  v  ltage  is  applied  to  the  circuit.   If  an  initial 
condition  is  prescribed  appropriately,  a  subharmonic  oscillation  of 
20  cps  may  be  started  in  the  circuit.   If  one  theoretically  used 
a  transformer  core  wh<  >se  characteristic  is  expressed  as 

3 
f  (v)   =  c,v  +  c  ,v 
1       3 

the  graph  of  Figure  2-2  is  i  btained.   However,  the  nonlinear 
characteristic  of  the  ordinary  transformer  core  is  not  truncated 
with  the  cubic  term  and  is  expressed  by 

3        5        7 

f(v)    =    C,V   +   C -.V     +   CCV    +   C,V    +    «ooo 

N   '  1  3  5  7 

1    the   regi   n    m  which    1/3-harmonic    oscillation    is   sustained 
appears    Ln  Figure   2-1. 

By  connecting  a  number  of  inductance  coils  in  series  and  ad- 
justing the  length  of  the  air  gap  which  is  interposed  in  each  core 
a    fairly  g     d    approximation   of   the   characteristic 

3 


f(v)      -      c    v     +      c    v 


is   obt ained. 
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Figure  2-1.   Region  in  which  1/3-harmonic  oscillation  is 
sustained.   Magnetization  curve  by  the 

equation 

3 
f(v)=c1v+c„v  . 
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ire  2-2.   Region  in  which  1/3-hannonic  oscillation 
is  sustained.   Magnetization  curve  by 
equation  3    5    7 

f  (v)=C.V+C3V  +C5V  +C.-V  +.... 
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When  two  cores  are  used2  one  with  an  air  gap  and  the  other 
without s  the  experimental  verificati  n  to  the  analytical  analysis 
is  quite  satisfactory „ 

2,1.2  Experimental  investigate    t  subharmonics  f  order  1/2, 
Investigations  of  systems  described  by  the  equations 

4k,  *  28    t£><^  =  Beet  ZU>t*   B, 

dl  4?  (2-9) 


*/t  ^  (2.10) 

can  be   done   using   an  analog  computer.      Figure    2-3   shows    the  block 
diagram  of  an  analog   computer    set   up  of  Equation   (2=10).      To  obtain 
a   set   up  of  Equ   ti        (2-9)    one   reeds   to  replace   the   servomultiplier 
by   an   ordinary  multiplier  which  gives   a  cubic  nonlinearity. 

If   the    initial  conditions   of  v(0)    and  v(0)    are   chosen  appropri- 
ately one  may   start   1/2-harmonic   oscillations,      Slowly  varying   the 
vi  ies      f   B   and  B     produces   the   regions    in  which   1/2-harmonic s   are 
sustained  and  Figures    2-4  and   2-5    show  the   experimental   results    for 
Equations    (2-9)    and   (2-10)    respectively. 

Figure   2-6    is    a  graph      l   the   region   in  which   1/2-harmonic 
oscillatio   s   are    sustained  using   theoretical  calculations. 
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Figure  2-3.   Computer  block  diagram  for  Equation  (2-10) 
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Figure  2-4.   Region  in  which  the  1/2-harmonic  is  sustained; 
nonlinearity  by  cubic  function  (analog -computer 
analysis) . 
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Figure    2-6.      Region   in  which  the    1/2-harmonic    is   sustained 
(calculated) . 
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One-half  harmonics  are   also  witnessed   in   the   electric  circuit 
of  Figure    2-7   due   to  the  presence   of  a   saturable -core    inductor  L 


-V 
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*» 


<C 


1  (- 

-*•      C 


r\j     )    E  Sin  Z  wt 


Figure    2-7.        Oscillatory  circuit   containing 
reactor  with  direct   current 
superposed , 


under  the    impression  of   an  alternating  voltage      cScftiwt 
The    secondary  winding    is  provided  on  the   core   so  that   asymmetry   is 
afforded   to  the  nonlinear  characteristic  by   forcing   a  constant  DC 
current    flow  through   it. 

Using   the   notation  of  Figure    2-7  we  have   the  equation 


d  t 

Re,  -  Vc  [  tc  A1 


(2-11) 


where   n   is   the   -rumfeer    of   turns   in   the   primary  winding   and   0   is   the 
magnetic    flux   in  the   core. 
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The   nan  dimensional  variables   u,    u    ,    and  v   are    introduced   in 

.> 

place       f    is    i    3    and    0   by    che    re  1     w^-i 

(2-12) 

where   "I,,      and       <^M        .re    appropriate   b-*se   quantities   of   the   current 
and   flux0      Neglecting  hysteresis   the    saturation  curve  has   the    form 

3  5  7 

u     +      u       =     c   v     -!-     c  ,v       +     c-V       +     c    v       +   .0  0  (2-13) 

The  base  quantities  3^  -j;id  pn    may  be  fixed  by  the  relations 

By  eliminating  i  and  i   in  equations  (2-11)  one  obtains 
where   t«  Cut-  %    *«.*"'   */x       ^ ■*   Vw«C 


A 


I n      rier    that    the   experiment   more   closely  agree  with   the   analysis 
a  composite    reactor   as  described   in  section   201„1    is   used?    the 
saturati  n  characteristics   of  which   ire   described   by 

Thus    the    differe -,ci-l  equ^ciun   of    the   electric    circuit   becomes 


dV    4    ^    chr  _^  ^^  f  ^    „    ft    ^cos  2  tot*  ■♦ 


c 

i  C 


d 
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Figure  2-8  shows  the  experimental  data  obtained  by  varying 
the  DC  current  in  the  secondary  winding  and  the  applied  AG  voltage 
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Figure   2-8  „        Region  in  which  the   1/2-harmonic 
is   s us t a i ne d   (,  e  x pe r i  me u t a  1 )  , 
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2.2   Subharmonics  generated  in  a  closed  loop  system  with  thyristors 
amp  1  i  f  ie  r  s  , 
It  has  been  found  that  the  problem  of  subharmonics  is  not 
limited  merely  to  those  cases  described  in  the  foregoing  discussion. 
Dr„  F.  Fallside  and  Dr,  A,  R.  Farmer  have  investigated  the  occurrence 

i.  f  subharmonics  in  systems  which  employ  thyristors  to  create  a 

11 

power  amplifier,. 

Power  amplifiers  that  use  thyristors,  thyratrons ,  magnetic 
amplifiers  and  mercury-arc  c  iverters  are  well  known.   Like  other 
types  of  controlled-rectif ier  amplifiers s  the  thyristor  amplifier  is 
a  discontinuous  element  since  the  input  signal  to  the  amplifier  con- 
trols the  output  at  discrete  instants  of  time  only.   The  output  wave - 
fi  rm  is  characterised  by  step  changes  in  voltage  and  hence  is  rich 
in  ripple. 
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Figure  2-9.   General  c  ntro]  1   p, 


When  used  i  i  a  c  itrol  loop  such  ^  Figure  2-9,  the  thyristor 

be  most  simply  treated  as  a  c  ntinuous  element,  represented  by 
( 

T 

its    Incremental  DC   gain  between   its    steady-state   input   and    its  mean 
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outputs   and   the    ripple  voltage   due   to  the   output  wave   can  be   dis- 
regarded,,    This   is   satisfactory   for   low  bandwidth   systems  where   the 
maximum  component    frequency  of  amplifier   input    signals    is  much   less 
than   the   amplifier   sampling   frequency „     Alsos    the    large   time  constants 
i  ic    t'f'   rated    nt      the   linear    filter   is   sufficient   to  ensure    that 
only  a  negligible  amount   of   ripple    is   fed  back. 
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Figure    2-10, 


Waveforms   of   2~phase   power   supply, 

stable   operat  1.    i . 
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In  wide   b  a  ;d   ^ppli      I  Li    ls    such    as    curre    t-control    loops    the 
simple    approach    is    not   valid,       [n    such   systems  enough  gain  must   be 
used   si'    that   appreciable    ripp;e    is    fed  b^ck   into   the   amplifier    input 
a  id      u    ler    ze.x\    In  c     iditii  is,   a   self   oscillation  of   the    loop, 
known   as    ripple    instability   can   result.      The    frequency   of   the 
oscillation  occurs   at    an   integral    submuLtiple   of   the  basic   ripple 
frequency      i,e,,    a.    subhi,rmc  uc  .     -J    m  ■  lifests    itself    in   the    amplifier 
as   a   periodic   variation    in   the    firiag    pattern. 

The    wavef    rm  of  Figure    2-10   shows   the  behavior   of  the    typical 
stable    operation  of   a   2-phase   regulated  power   supply. 

Figure    2-11    Shows    the    ripple    instability    introduced  when   the 
loop  gain   is    increased  by    the    preamplifier. 
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Ripple  instability  cannot;  be  predicted  by  a  simple  linearized 
approach  to  system  design  -  d  a  more  accurate  analysis  must  be 
used  . 

2  . 2 . 1   Amplifier  characteristics 

A  thyristor  amplifier  may  be  generally  represented  by  the 
m-phase  arrangement  of  cells  shown  in  Figure  2-12.  The  output 
voltage  waveform  of  the  amplifier  for  a  constant  input  y  is  assumed 
t   be  as  shown  i   Figure  2-13  for  a  resistive  or  diode  clamped  load. 

There  are  many  types   f  firing  circuits  for  thyristor  amplifiers 
and  the  sinewave  reference  type  is  c  nsidered  here.  This  type  of 
firing  has  an  advantage  in  that  it  produces  a  linear  DC  gain  character- 
istic for  the  relation  between  the  amplifier  input  y  and  the  mean 

output  V'  ltage  v  .   In  this  particular  circuit,  the  input  to  the 

1 1 

firing  circuit   f  each  phase  consists  of  the  amplifier  input  plus  a 
phased  sinusoidal  reference  signal  of  constant  amplitude  at  the 
excitation  frequency.   A  gate  pulse  is  passed  to  the  thyristor  when 
this  total  input  passes  through  zero  going  pusitLve. 

The  Fallside  and  Farmer  papei  precedes  to  analyze  this  parti- 
cular amplifier  by  two  means9  an  impulse  method  and  a  Master 
Inscribing  F    L  i   metr!  i. 

2.2.2.  Analysis  by  impluse  meth  1 

Si   e  the  subharmonic  of  i  rder  1/2  is  the  simplest  oscillation 
t   analyze  by  tr  j  impulse  metr  d  s  that  it  only  will  be  considered 
here. 

The  u  b      l-amplifier  output  waveform  characteristic  of 
1/2  i  rder  subharm  i «.  instability  is  shown  in  Figure  2-14a.  This 
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Figure    2-12.      Generalized  arrangement   of  m-phase   amplifier 
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Figure    2-13.        Output  waveforms   of  m-phase   amplifier, 
resistive   or  diode -clamped   load. 
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a.      Unbalanced  waveform 
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Figure  2-14.  Unbalanced  operation  of  a  thyristor  amplifier 
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waveform  is  considered  to  have  two  parts;  the  balanced  component  of 
Figure  2-14b  and  the  palse  '  .•  i   >f  Fig. re  1  -  14c ,  alternating  at  the 
subh-.rrn  nic  frequency  if  mil  /2e   The  pulse  train  which  is  representa- 
tive of  the  imbalance  has  no  net  DC  output  and  injects  a  signal  of 
frequency  mXl  /2  Into  the  loop.   A.  portion  of  this  signal,  attenuated 
and  smoothed  b»  tre  filter,  is  fed  back  into  the  amplifier  input  and 
the  oscillation  m - y  thus  be  maintained.   F  r  any  appreciable  oscil- 
1  - 1  i  n  -b  it  the  mean  firing  angle  <*■     ,  the  positive  and  negative 
pulses  are  different  b  th  in  amplitude  and  duration.  However 9  for 
dimini shingly  smaller  amounts  of  imbalance,  the  pulse  train  becomes 
in  the  limit;,  a  symmetrical  impulse  train  of  Figure  2-14c,  and 
represents  the  conditions  at  the  onset  of  instability.   By  examining 
the  system  equati  ns  under  these  c  md it ions  it  is  possible  to  formu- 
late a  general  c  nditi  n  for  the  loop  to  sustain  a  subharmonic  of 
infinitesimal  amplitude. 

Referring  tc  the  linear  c  ntr<  L  1   p  of  Fig.  2-9,  assuming 
stables  balanced  sperating  c  mditions,  the  amplifier  output  can  be 
written  as  the  Fourier  series 

The  input  t   the  thyr is t  r  amp L if ie r  is 

'        '  l  }  (2-15) 

and  the  equation   f  the  firing  circuits  at  any  firing  instant  is 

>  (2-16) 
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whe  re 


Jit=«  +  1^-  (k*o(i,v  .  . ,-)  (2i7) 


and 

p    -  1 


f   s  ili.  +  ff  . 


(2-18) 
Substituting  f  r  y  Equati  n  (2-16)  becomes 

K^a;-   K^o)^(o)u-  +  2cos(^-^-r)=o        (2_19) 

Considering  steady  state  operati  n  where  v   is  constant,,  the  operator 
D  may  be  replaced  by  j  U)    ,  after  substituting  for  £  ,  SL~k      ,  and  P 
Equati  a   (2-19)  becomes 

n-.  i 

(2-20) 

-  j£  COS  ex   =  O 

The   operator     %(&)     and       </^0)       have   been   replaced  by    ^(j^)   and 
and   7^(j'u))  xs    the    amplitude   and  6?^      the   ph-*se   of  the  expression  at  a 
frequency    60-  nmil  „ 

If     ne   now  assumed  an   infinitesimal  amplitude   of   subharmonic 

sell lat ion,    the    system  equati     *s  must   be  modified  by  the   pulse   train 

£>\j-    ,    assuming    the   pulses  hive   a  duration      Za°<> 

Alt:    £   ^    *»%  »"«  CM  U  (At  -%  -  tgi) 

£.'/»  (2-21) 


where 


J?z     H-  '/2  (iflsl,?,"  ',  o°) 
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The    fir  leg    Lnst       ts     low  become     OCX-  (^-^Qt  )■+  -J^T"    ;    thu 
the    firing   circuit  equations   be  ernes 

-l\2,    *<        ^       ■    stn~  sen*  Cos  )*»■*(    ~ 


(2-22) 


-     A*  +  Hi  )  +   ^  ,0j-  J  CCS  (oU  ^  oc)  a.    O  . 

Assuming  /V-<  -•  0  ,  the  terms  which  represent  balanced  operation 
may  be  removed  by  subtracting  Eq0  (2-20)  „   The  conditions  by  which 
subharmonics  may  exist  is  seen  ^s 


1 1 


-   Zl   2n  Cn(°<)  n  m  &*  Cos  (n  iyic(  -*  <^.  *  <:-n)/ 

r>-  t 

Dividing   through  by    l^^i^'A      ,    and   introducing   the  DC   gain 

,.         iry*£  7T 

k  .  -    -_  St  n  "TT «      t  h  is   reduces    t 


(2-23) 
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This  may  be  written  more  simply  as 

OO  oG 


Hx  "=l  (2-24) 

where  Q      =     Cn(*)t\rA_ 


X  ICA  CO%(oi-V^) 


a  i  j 


y  =  nm».-'  s^n  . 


In  summary,  the  critical  Loop  gain  for  ripple  instability  at 
a  subharmonic  of   rder  1/2  can  be  calculated  from  the  impulse- 
analysis  criteri  n  of  Eq.  (2-24)  for  any  given  thyristor  amplifier, 
in  any  linear  control  system  to  any  degree  cf  accuracy,  by  taking 
sufficient  terms  in  the  summation.   In  Eq,  (2-24),  C  ~  and  0_   depend 
the  amplifier  used  ->nd  its  mean  firing  angle  or  output  voltage; 
the  terms  x-   ,  •£?    ,  6j»   and  S      define  the  remainder  of  the  control 
1  iop  thr  ugh  ^(juj)    Yt C  j  Oo )  » 

2.2,3  Describing  functi  i  for  infinitesimal  subharmonic  input  signals 
of  t/rder  1/2. 

'The  Fallside  and  Farmer  paper  also  presents  a  describing 
function  approach  to  the  analysis  c £  the  conditions  for  sustentation 
of  subharmonics  of  urder  1/2  considering  the  subharmonics  appear  as 
infinitesimal  amplitudes.  Again  f  r  a  constant  input  the  firing  cir- 
cuit  equation  is 

tlcoa  (fL-t  -Vm  -e)-o 


* 


(2-16) 


id   at    t  he    firing    inst  ant  s 

nt -«**£*     (k-v-v,-) 
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If   -i   si-.as    id    1    sig  ,.L   of    infinitesimal  amplitude   LSu     ,    frequency 
rr>Xly^  i   some   arbitrary   phase  u     is   -died   to  y,    the   firing 

equations   become 

m  +    IS.  cos  (  —>  +  ^)  +  2  Co5(jftt  -  %      |   )  =  O 

and 

jit   =  ^i  t  &*  +    T?r 

After    subtracting   Eqn„    (2-16),    the    incremental   firing   angle 
•i-id    the    i    put    subharm   nic    signal  amplitude   Aij     are    related  by 

=  -  A-Z^>(**  \vr-£-f).  (2-25) 


The    i^Jici    n  of  a   subharm    lie    s  i  fc.  -  <  -.  1    t  i   the    input   of   the   am- 
plifier   in   turn    produces    an  alternating    impulse   train  at   the   amplifier 

jtputo    iFi^„    2- Lid)        This    impulse    train  has    a  component   of    sub- 
harmonic    frequency  given  by 

(j-y    =    -—  A^  Sen  /W|5*n«i  Cos  l    L    t-_ .  -TTkJ 

Substituting    from  Eq„    (2-25)    f    r    Au     and  assuming     A<?<  <  <  ^  j 
Uy      =    £  ICA  COS     ^  CO 5    (  %-  *    +  /a   -    00 


=    ZK<  CO)     V  Co 5    ^ 
where  (2-26) 


Thus ,    f    r    i    pot    sign    is      E    i--.fi  Litesima]    amplitude   at   subharmonic 
frequency    r!3^        ,    the    amplifier  has   a  g^in  of  ZK/tCosfP  and 
intr  duces    a   phase    shift    ^  . 
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The  amplifier  cam  then  be  characterized  by  the  describing 


function  § 


NOf 


n-l     ^a 


eoi(^yu) 


2KxCo^\^ 


(2-27) 


This  is  shown  in  Fig.  2-15 9  and  is  a  circle  of  radius  Ka   ,  centered 
on  (  K»  \l  0  )«  f'he  normalized  version  [Sly,*  ~T~/   LS  a  circle  of  unit 
radius  centered  on  (l9jO). 

This  describing  function  may  be  used  to  test  for  ripple  in- 
stability at  the  subharmonic  frequency  ~r—   and  the  procedure  can 
be  seen  in  Fig.  2-16.   First  the  open- loop  frequency  locus  l\K$  **,(}  <*mi(}u>) 
is  plotted;  then  the  normalized  describing  function  Nvj  "7")    Is 
superimposed j  as  a  unit  circle,  centered  upon  the  point  on  the 
frequency  1.  cus  corresponding  to  CO-  ~^~   .   Using  conventional 
assumptions,  the  system  will  be  unstable  and  oscillate  at  the  sub- 
harmonic  frequency   — r—   Lf  the  circle  encloses  the  point  (-l,jO); 
the  critical  condition  being 

2KKA   £  1/2  cos  0  l/2   ■»  -1    =0  (2-28) 


Referring  to  Eqn.  (2-24),  this  is  the  condition  which  was  obtained 
from  merely  the  first  term  of  the  exact  impulse  criterion. 
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Figure  2-15. 


Describing  function  for  thyristor 
amplifiers  with  subharmonic  input 
signals  of  order  1/2. 
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Figure    2-16.      Use   of  describing    function  to  test   for 
ripple    instability. 
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3,  Description  of  system 

In  chapter  2,  the  SCR  amplifier  that  was  discussed  was  one  that 
employed  a  phase  firing  scheme  of  control.   The  thyristors  were  fired 
individually  when  the  current  ir  r  jgh  each  one  attained  a  certain 
pi  se  in  relation  to  a  reference  sinusoid.  The  thyristor  amplifier 
that  is  the  concern  of  this  paper  employs  a  different  method  of 
operctt  ion. 

The  specific,  system  under  i  ivestig-ili  n  is  a  polyphase,  thy- 
ristor, AC  to  DC  rectifier 9  operating  under  forced  limit  cycling 
conditions,   The  block  diagram  for  the  model  of  the  system  appears 
in  Figure  3-1. 
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Figure  3-1.   Bl  ck  iiagram  for  Regulated  Thyristor  Power 
Supply  - 


This  block  diagram  may  be  simplified  into  the  system  of  Figure 
3-2  where  G-  is  Included  to  take  into  account  any  gain  that  may  be 

L 

s   ciated  with  the  comparators  firing  circuits  and  compensation  that 
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may  be   -?,dded   to  the    system  and   is   completely    linear,      G      is   the 
linear   transfer   function   of   the    filter   and   loadc      All   of   the   non- 
linearities    of   the    finig   circuit    a  nd   bridge    r,  re    lumped    into   the 
single   describing    functi   n  N„      The         ture     if    this  describing   function 
is  developed    in  Leszcynski    s    thesis   not   considering    ripple    instability 


10 


F  1  g  3 ;  e    3 - 2 o         S imp li.f  l si  b  1  oc k  d  i  ? g x  -?m < 


In   the    pr^se    firing   scheme,    for  polyphase    rectifications,    It 
is    required   that    there   be    separate    firing  circuits    for  each    individual 
thyrist    r   as   well  as    phasing   and   timing   circuits   to  ensure   proper 
sequence    of    firing.      Hence,   a   three  -phase  ?    full  wave    rectifier  would 
require    six   separate    firing   circuits., 

So   ^s    '      simplify   the   c.   ttrol    phi]  isophy   and  to  curtail  the 
expense      f    six   separate    firi-ig    arc  iits9    the   method  considered  here 
is    c      n-.tr    i    the   entire    p    Lypt    se    bridge    as    a   unit   with    one    firing 
circuit,      1      this  way    the   g->res      f   -.  11    the    thyristors    in    the  bridge 
-^re   c    tinected    in    series    and   are   pulseH    simultaneously  with   one    pulse  „ 
(A   simplified   circuit    diagram  can  be    seen  in  Figure    3-3„)     However9 
only  those   thyrist    rs    ah    se   current   has   a  phase  between  0  and   180 
degrees  will    fire    immediately   as   the   pulse   gees  ON.     As   the   phase 
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of  the  current  Lm  the  remaining  thyristors  traverse  through  0  degrees 
they  too  will  fire  as  lung  as  the  ON  pulse  is  maintained  at  the 
thyristor  gates.   Figure  3-4  illustrates  the  sequence  as  the  six 
thyristors  are  pulsed  ON.  As  a  result9  when  the  firing  signal  is 
ON,  the  entire  bridge  is  on  and  acts  as  a  conventional  diode  bridge. 

A  slightly  different  aspect  is  presented  when  the  bridge  is 
pulsed  OFF,,  or  the  bia,s  is  rem  ved  from  the  thyristor  gates.  When 
the  OFF  pulse  is  received,  a  time  delay  occurs  before  the  bridge  is 
completely  cut  off.   This  is  due  to  the  fact  that  the  thyristor, 
like  its  gas  tube  counterpart;,  the  thyratron3  is  not  cut  off  until 
its  shut  off  conditions  have  been  met,  i.ec,  the  bias  has  been  re- 
moved and  the  collector  voltage  goes  negative  with  respect  to  the 
emitter.   Figure  3-5  shows  the  sequence  as  the  six  thyristors  are 
pulsed  OFF.   The  time  delay  C   occurs  between  the  time  the  OFF 
signal  is  received  and  the  last  thyristor  cuts  off. 

A  summary  of  the  system  operation  is  as  follows.   Assume  that 
the  desired  output  is  a  certain  voltage  and  the  output  voltage  is 
below  this  reference  level.   The  rectifier  bridge  is  then  supplied 
with  an  ON  s ig  al  and  the  bridge  output  is  a  three-phased,  full-wave, 
rectified  voltage  which  is  fed  to  the  filter.   The  filter  then  at- 
tenuates the  harmonics   f  this  voltage  and  provides  a  relatively 
smooth  DC  voltage  at  the  output.   When  the  DC  output  rises  above  the 
reference  level  the  bridge  is  pulsed  off.   Because  the  shut  off  con- 
ditions of  each  thyrist  r  most  be  met,  a  time  delay  occurs  before 
the  bridge  turns  off.   I  i  the  OFF  state  the  input  voltage  to  the 
fitter  is  essentially  zero  since  the  reverse  current  is  virtually 
blocked  i  ff  by  a  clamping  diode.   The  briige  thus  remains  in  its  OFF 
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Figure    3-3.        Three-phase,    full -wave,    forced   limit -cycle 
regulator. 
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Figure    3-4.        Firing   sequence   of   thyristor  bridge. 
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Figure    3-5.        Cut-off   sequence   of   thyristor  bridge. 
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state  until  the  output  voltage  drops  below  the  reference.,  whereupon 
the  cycle  repeats  itself. 

Leszcynski's  thesis  shows  that  the  output  of  the  power  supply 
is  a  DC  voltage  with  a  ripple  that  is  a  non- sinusoidal  limit  cycle. 
with  a  fundamental  frequency  different  from  the  AC  supply  voltage  or 
any  of  its  harmonics. 

The  input  to  the  filter  can  be  shown  as  a  series  of  pulses, 

one  of  which  appears  in  Figure  3-6  below, 
Input  to  filter  (volts) 
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Figure  3-6.    Thyristor  bridge  output  pulse 


Suppose  that  the  period  of  the  supply  voltage  is  T  seconds.   By 
examining  figure  3-6  it  can  be  seen  that 3  depending  upon  the  length 
of  the  firing  pulse,  there  is  a  period  of  T/63  the  time  that  exists 
between  one  thyristor's  firing  and  the  next  thyristor9 s  firing,  in 
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which  the  cut  off  of  the  firing  pulse  will  not  vary  the  elapsed  time 
of  the  output  pulse  of  the  thyristor  bridge.   It  is  also  seen  that 
the  time  delay   C   can  vary  from  T/3  to  T/2  seconds  in  length,  or 

T/3<  t   <      T/2 

It    is   this  variation   in     C     that  might   give    some    insight    into 
the   problem  of    ripple    instability   in  a   forced    limit-cycling   regulator 
and  will   be   examined  more   carefully   in   the   next   section. 

3.1        Preliminary   investigation   into   ripple    instability. 

It  was    seen   in  the    preceding    section  that   the   output   of   the 
thyristor  bridge    under   steady-state    operating   conditions    is   a   series 
of   pulses   that   take    the    form  of   the   pulse    shown   in   figure    3-6. 

If    one    assumes    that    the    ripple    of    the   DC   output   of    the    filter 
is   a    limit-cycle    that   closely  approximates   a    sinuooid,   with   a    funda- 
mental  frequency   f      and  a   fundamental   period  T    ,    then   the    length   of 

the    firing   pulse   would  be  I\/2  or   the    length   of   time   that  V      is    less 
°    r  1  °  o 

than  V      in  one    limit -cycle   period  T 

Figure    3 -7a    is   a   replica   of    figure    3-6   in  which    the    firing 

pulse    is  T./2   seconds    long   and   the   thyristor  bridge   pulse    length    is 
i 

t1    seconds.      Figure    3-7b    shows   the    firing   pulse    still   to  be   T/2 
seconds   but   displaced    in   time    from  the    first    firing   pulse   by  € 
seconds.      The    resulting   pulse    formed  by   the   thyristor  bridge    is   t 
seconds  which  can  be    seen  to  be    less   than  t..  , 

Suppose   one   used   this    sine   wave    limit   cycle    over  a    long   period 
of   time    to  generate   a    train  of   pulses    that  might  be   obtained   from 
the    thyristor  bridge    under    steady-state   operating  conditions.      Assume 
that   the   bridge    is   pulsed   on  when   the    sine  wave    is   greater   than   its 


45 


Input   to   filter 


/N        / 


/    X  '    \ 


1 i )L v  v 


^ i. 


Firing 
Pulse 


Tx/2 


*\ 


Figure   3-7a. 


Input  to  filter 


\  1 


V      /        \7     w       y    '   \ 

v  /  /  \  \     I 

'*     /  \     A     A     /  \       v 

; ¥ \l v \       \ 


F  ir  ing 
Pulse 


\, 


Figure    3 -7b. 


46 


N 
4-1 
U 

0) 

m 


o 

o 
I 


CO 
P-4 


CU 


CO 

> 

CO 

13 
ti 
O 
cj 
CD 

CO 


5-1 

£1 

o 

CU 

4-1 

CM 

j: 

T3 

i-i 

& 

<t 

CU 

ii 

CO 

.—1 

CU 

3 

r-l 

Ph 

O 

P*> 

O 

O 

r-l 

1 

o 

■P 

00 

•H 

1 

II 

CO 

hJ 

CU 

u 

00 

• 

■r-l 

CM 
O 


0> 

i-H 

o 
o 


.-i 
o 

• 

o 


CM 

o 


ON 

T-l 

o 
o 


47 


N 

u 

0) 
X! 

CO 
CO 


o 
o 


•H 

►4 


CM 

O 


m 

CM 

o 


48 


o 


w 

f> 

o 

w 

CO 

1 

O 

o 

00 

• 

c 

N 

•^ 

•U 

u 

J3 

o 

0) 

■U 

CM 

J3 

•H 

m 

& 

CM 

<D 

II 

CO 

<u 

3 

r-t 

PM 

o 

>•> 

O 
<-l 

o 

1 

- 

4J 

00 

•H 

1 

CO 

l-J 

a) 
60 

• 

•H 

0) 

Pn 

CTi 
CM 
O 


1^- 
CM 
O 


CM 

O 


49 


mean  value    a:    1    is   shut   off  when  the    sine  wave   goes  below   its  means 
s      as   ti     simulate   the    freqp.  icy      f   the   error  signal  becoming  greater 
or   less   than  zero.      Observed   ever   fifty  or   sixty  cycles   of   the 
frequency   f11    one   notes   that   the  pulse   train   from  the  bridge   t-jkes 
i    an  >hvious   pattern  as  can  be    seen   in  figure   3-8a-e, 

Observe   that   the   period   of   repetition  of   the   pattern  established 
by  this  means    is   some   times  many   times   greater  than  the    limit   cycle 
peri  od  T.    but    In  all  cases    is   an    integral  multiple   of   the    limit-cycle 
period.      If   one  were   to  write   a  Fourier  Series   representation  of 
this    pulse    train,,    as   complicated   as    it  might   be,    the    fundamental    fre- 
quency  in  each  case  would  be   a   subharmonic   of   the    sine  wave   that 
generated   if.      Hence s    it    can  be    seen  that   a  pulse   train  generated  by 
a    si nus    Ldal   signal  with    frequency   f-,    presents  energy  to  the    filter 
that   contains   subharmonics   of   the   generating   signal. 

This  s    of   course.,    is  not    irrevocable   proof   that    subharmonics   do 
exist    in   the   output  because    it    is   known   that   the    limit-cycle    is  not 
a    true   sinusoid  and    also    >ne   does  not   knew  what   farm  of   pulse   train 
is  generated  by  a   limit -cycle    in  which   ripple    instability  does  exist  , 
but    one  must    admit   to  a  definite   possibility  that   subharmonics   of 
the    limit-cycle    frequency  might  be   present   in  the   output. 

The  next   step  would  be    to   simulate   the    system  on  an   an* log   or 
digital   computers    if  not    the    actual   plant    itself,    and   in   fact    see 
if   subharmonics    ^.re   present    i      the    i  utput   and   under  what   conditions 
they  exist „      This    is    indeed  what    will  be   done    in  the   next   chapter. 
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4.   Synthesis  of  the  forced  limit -eye ling  regulator. 

The  system  under  study  was  successfully  synthesized  in 
Leszcynski's  thesis  by  a  digital  computer  program  utilizing  Fortran 
IV  language  and  the  Naval  Postgraduate  School  IBM  360  Computer  Facility 
In  the  program  the  firing  circuit  and  SCR  bridge  was  simulated  by 
logic  statement  manipulation  and  the  output  of  the  linear  filter  was 
obtained  by  calling  a  library  subroutine  INTEG  2  which  gives  a  Runge- 
Kutta  solution  to  a  group  of  differential  equations,  in  this  case 
the  state  variable  equations  developed  from  the  dynamics  of  the  linear 
filter  (see  Appendix  I) .   The  linear  filter  is  assumed  to  be  as  shown 
in  Figure  4-1,   The  transfer  function  for  the  filter  and  load  is 
developed  as  follows; 


(J-    =  (TL  +  0c,         t  L  =  v.  c  +■  <-  R 


rlr.  °  cit  R 


-  L-r^t 


—  -I  I  in    <-)        dt     C     RC 


where    tC  —  l-"^ 


Choosing  the  state  variables  as 

a  signal  flew  graph  of  the  linear  system  is  shown  in  Figure  4-2. 

Using  Mason's  gain  rule  the  transfer  function  in  Laplace 
notation  is: 


T(S)   = 
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Figure  4-1.   Filter  and  load. 
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Figure  4-2.    Signal  flow  graph  for  linear  system. 
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The    feedback  path  between  the   error   signal  and   the    firing 
circuit    is  known   in  the   actual   system  to  contain  a  dead  zone „      Figure 
4-3   is   the   dead   zone   that  was   assumed   in   simulation  of   this    system, 
with     CX    being   the  width  of   the  dead   zone.      In  the   program  a  gain  C(6) 
is   introduce  1    Lnt  )   the    feedback  path    t      pr   vile    amplif icationrto  the 
error   signal   prior   to   the   dead   z  •  Le       id   has   the    ability   of  defeating 
the   effects   of   the   dead   zone    if   used   properly  as  will  be    seen   in  the 
experimental  results. 

In  conducting   the   experiment   the    supply  voltage  was   considered 
to  be    a   three-phase   60  cycle  AC    source  with    a   peak  amplitude   of   65 
volts.      The   value   of   the    inductance   and  capacitance   of   the    linear 
filter  were   kept   constant   throughout   the   experiment   at   0e0005  henries 
and  0,1   farads   respectively* 

Thus s    the   experimental  data  was   collected   using   the    following 
factors   as  variables: 

C(4)  -  The   value   of   the    load   resistance    in  ohms 

C(5)  -  The   value   of   the   reference  voltage 

C(6)  -  The   amount  of  error   signal  amplification 

C(7)  -  The  width   of  the   dead   zone    ( 0<  ) 

Using   a   step   size   of   0.00002222    sec    nds    it  was    found   that    in  most 
cases  by   setting   the    initial  conditions   close    to  the  desired  reference 
voltage    the   transient   died  out  within  the    first   0,1   second   and   the 
effect    of   ripple    instability  could  be    successfully  observed   from  0,1 
tc    0,2   seconds.      Thus.,    the   graphs   that    ^re    shown  begin  at   0,097687 
seconds   and   continue   to  0.19702   seconds.      The   effects   of  manipulating 
the    sever -il  variables   can  well  be    seen  during   this   time    frame    in 
most    cases. 
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Figure  4-3.   The  dead  zone  characteristics  of  the  feedback 
system. 
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4.1  Experimental  results 

The  first  c  ncern  was  to  find  mt  if  the  output  of  the  filter 
actually  contained  some  subharmonic  ripple.,  s   taking  time  intervals 
of  Ool  seconds  the  simulation  program  w-s  run  until  the  transient 
had  obviously  died  away  and  the  steady  state  output  was  examined., 
Fig,  4~4b,  sh  iws  the  steady  state  output  of  the  regulator  with  the 
load  equal  to  1,6  >hms  and  the  reference  voltage  set  at  60.0  volts, 
(The  dead  zone   in  each  case  is  set  at  0„0  volts  unless  otherwise 
noted  on  the  graphs „) 

The  load  resistance  was  subsequently  doubled  and  then  halved  to 
see  the  respective  results  on  the  output.  Those  results  can  be  seen 
in  Fig,  4-4a  and  c.   Fig.  4-5  shews  the  input  pulse  train  for  the 
corr e s pond ing  t  ime  pe r i od 0 

In  Fig.  4-6a  the  reference  voltage  with  the  load  set  at  3,2 
ohms  was  dropped  to  58,0  vults,   Fig,  4-6b  presents  the  results 
when  the  err.  r  signal  vis  doubled  but  zero  dead  zone  was  considered. 

Keeping  the  I  ad  resistance  at  0,8  ohms,  the  reference  voltage 
was  dropped  in  steps  from  65,0  volts  to  48,0  volts.   These  results 
appear  Ln  Figs,  4-8s  4-93  and  4-l0„ 

Using  the  values   f  R=0,8  ohms  and  V  f  =60„0  volts  the  system 
was  examined  assuming  the  dead  zone  to  be  0,5  volts.  These  results 

r re  compared  t  >  the  results  obtained  by  setting  V   f=60,5  volts  and 
the  dead  z   e  set  at  0,0,  Fig,  4-Tl  shows  this  comparison. 

Figure  4-13  is  a  c  mp  ris  n  between  the  results  btained  by 

setting  V   =62,0  volts  and  the  dead  z  ne  at  0,0  and  setting  V   .=60,0 

°     r  e  i  re  r 

v  Its  and  lead  zone  equal  tc  2,0  volts.  The  err  r  signal  was  then 
amplified  by  a  factor  of  L.5  and  compared  to  the  proceeding  experiments 
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Finally  the    system  was  examined  when  the   pulse   train  was 
amplified  by  a   factor  of   2.0  before  going   into  the    filter  and   the 
results  were   compared   to   setting   the   source  voltage   to  double    its 
original  value.     These  experimental   results  appear  in  Fig.   4-14, 

401       Summary        3    malysis   of    the   experimental  data. 

The   transfer   function   for   the    linear   filter  as   stated  previously 


is; 


1/LC 


T(]oJ    )   =     — 

(j  CJ)"     +      (jOJ  )/RC      +   1 

En  the   Bode    form  this  becomes 9 

1 


T(jO)  )      = 


',    (j  U>)         +      (J  <*>)L/R     +      1 


In  dec  i.bels 


T(j«*>d^     =   -20    I    g 


LC(j  U)  )        +     L/R(jCJ  )     +   1 


.    mparing   this   to  the  normal  form  »,  f  the   Bode  plot   equation 


T(j  U)  ) 

db 


=   -20   1   2 


(j  °°  /oun  )2  +  2  ^  /ujn  (J«o)  +  1 


it  can  be   seen  that 


<<->, 


V 


LC 


2  *    /mj         ■   L/R 


i  - 


1./2R 


s  i  thai  f  r  the  values  of  L  and  C  used  in  this  experiment 


H>  = 


141  rad/sec 


22o4  hertz 


/ 


( ,0005)(.l) 
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.  id 

1/2R 


C      -      1/2R          A2225          -          003535534 
"V     .1  "  R 


hence,    for 


R=3.2  C    =0,0110485 

R  =  l,6  £    =0.0220970 

R=0,8  £     =0.0441940 

Using   these    three   v-iues   of    load   resistance    the   B>de   diagrams 
for   the    linear   filter  response   c->a  be    seen   in  Figure  4-16. 

Table   4-1    is   a   tabular    summary   of    the    results   of    the   experiment 
in  which   the    reference   voltage   was   kept   constant   at   60,0  volts   and 
the   value   of   the    load   resistance  was  varied.      In  Table   4-1, 

R     =     Value   of   the    lead   resistance    in  ohms, 

f  =     Frequency    of    the    fundamental    limit-cycle    in 

he  r  t  z 

form      ~     Frequency   of   the   obvious    subharmonic    in 
SUB  . 

hertz 

Vp.   =  Pe^k  to  peak  amplitude  of  the  larger  spike 
i-i  the   utput. 

Vp.,  =  Peak  t   pe-.k  amplitude  of  the  smallest  spike 
in  the   utput . 


V   =  Damping  c  efficient 
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The  values  of  the  subhdrm  jic  frequencies  are  noted  for  each 
case  on  the  Bode  diagrams  of  Figure  4-16. 

Table  4-2  is  a  tabulated  summary  of  the  experiment  in  which  the 
load  resistance  was  held  constant  at  0,8  ohms  and  the  reference  volt- 
age was  decreased. 


R      =      C 

.8  Ohms 

fSUB 

V*l 

V?2 

V      , 
re  f 

^LC_ 

65,0 

27.56 

0 

10.004 

64,0 

32.91 

0 

8.001 

62.0 

35/79 

17.89 

6.927 

6.844 

60.0 

45.06 

22.53 

5,904 

4.192 

58.0 

51.79 

25.89 

5,003 

4.232 

56.0 

60.13 

15.03" 

4.997 

----- 

54.0 

51.79 

36.36 

3.301 

3.240 

52.0'wC 

60.93 

30.47 

4.410 

4.709 

50.0 

90.25 

45.13 

2.592 

2.551 

48.0 

99.99 

? 

■  - 

----- 

1/4-c 

>rder 

subharmonic 

Data 

does 

not    f 

1L 

w   p-* 

Item. 

Table   4-2. 

Because    the   physical   system   is  known  to  contain  a  dead   zone    in 
the    feedback  path,,   Table   4-3    is    a   comparison  of   the    system's   response 
when  a  dead    zone    is    assumed  and   the    system's    response  when  no  dead 
Ls  ci  ms  i  lered. 
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45.06 
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22.53 

5.904 
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5c   Conclusions 

The  main  conclusion  that  c  m   be  drawn  by  examination  of  the 
experimental  data  Is  that  under  certain  conditions  the  forced  limit  - 
cycling  regulator  discussed  will  support  subharmonic  ripple  instability, 
The  other  object  of  the  investigations,  that  of  determining  the  exact 
causes  for  the  subharmonic  sustentation,  was  not  clearly  achieved 
and  the  c<  aclusions  that  might  be  drawn  toward  this  end  are  somewhat 
speculative  and  will  appear  in  the  next  section  in  a  discussion  of 
the  recommendations  for  further  study  in  this  area. 
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6.        Recommended  areas   for   further    study 

One   tentative   conclusion  that  might  possibly  be   reached   in 
reviewing   the  experimental  results   is   that   the   subharmonic   that   is 
sustained    in  the   system  is  dependent   upon  the  numerics   of  the  com- 
ponents  of   the    filter   that    is  employed.      The    filter   in  this   case  was 
b    th   simple   and  a   little   unrealistic   but   it  did  bring   the    interesting 
observation  that   the    subharmonic   that  was   the  most   evident    in   the   out- 
put was  une   that  -was   quite   near   the   resonant   frequency  of   the    linear 
filter.      Examining   the    results   presented   in  Figure   4-4  one   can  see 
that   althdjtfgh   the   damping   of   the    linear   filter    is    increased  as   the 
i    ad   resistance    is    lowered,    the   amplitude   of   the    subharmonic    is    in- 
creased.     This  might  be  explained  by   the   Bode   diagram  of  Figure   4-16 
where    it   can  be    seen  that   as    the    load   is   reduced  the    limit-cycle 
frequency    is   also   reduced  and    its   1/2-order   subharmonic    is  moved 
closer  to  the   resonant    frequency  of  the    filter. 

In  Figure   4 -9b,,  where    the   resistance    is   0.8  ohms   and   the 
reference   voltage    is   56.0  volts,    the    1/4-order   subharmonic   of  the 
limit-cycle    is  more   evident   than  the    1/2-order   subharmonic   as  was 
the   case  with   the    rest   of   the    results.      This   also  could  be   due    to 
the   fact    that    this   order  subharmonic  was  closer  to  the   resonant  peak 
than  was  the   1/2-order  subharmonic. 

Based  on  the    above   philosophy,    a   recommendation   for   further 
study  would  be   to  refine   and  design  the    filter  to  produce   a  desired 
limit -cycle  whose    subharmonics   are  well  attenuated  by  the    linear 
filter. 

It  can  be   noted   in  Table   4-3   that  as   the   dead   zone    in  the 
feedback  path  was    increased   the    limit -cycle    frequency   decreased   and 
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its  1/2-order  subharmonic  moved  away  from  the  filter's  resonant  peak. 
Any  further  work  should  also  include  a  m  re  detailed  study  of  the 
effects  of  the  dead  zone  on  the  system. 

Another  area  for  study  that  was  touched  on  briefly  but  in  no 
extensive  detail  was  the  matter  of  varying  the  supply  voltage  and 
studying  its  effect  on  the  ripple  instability „ 

In  the  realm  of  mathematical  analysis  several  of  the  references 
noted  in  the  bibliography  suggest  that  a  possible  method  of  predicting 
the  perf  rmance  of  the  system  would  be  to  employ  a  Dual  Input 
Describing  Function  or  even  a  Multiple  Input  Describing  Function, 
Surely^  any  further  study  should  give  some  consideration  to  these 
analytical  methods  in  attempting  to  understand  ripple  instability  in 
this  system. 
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APPENDIX  I 

Computer  program  for  simulation  of  system  with  numerical  data 
for  the  conditions  when  the  reference  voltage  is  held  constant  at 
60.0  volts  and  the  load  resistance  is  3.2  ohms,  1.6  ohms  and  0.8 
ohms  . 
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£***** ***** ***** ***** ***** 

c 
c 

C  C<1»=MAXIMUM  VALUE  OF  INPUT  SINE  WAVE 

C  C(2)=VALUC  OF  FILTER  INDUCTANCE  IN  HENRIES 

C  C(3)=VALUE  OF  FILTER  CAPICITANCE  IN  FARADS 

C  C(4)=VALUF  OF  LOAD  RESISTANCE  IN  OHMS 

C  C(5)=REFE«ENCE  VOLTAGE 

C  C(6)=G£IN  OF  AMPLIFIER  BEFORE  DEAD  ZONE 

C  C(7)=MAGNITUDE  OF  DEAD  ZONE 

C  C(8)=GAIN  OF  AMPLIFIER  AFTER  DEAD  ZONE 

C 

C 

C***** ***** ***** ***** ***** 

c 
c 

DIMENSION  X<30) , XDOT< 30 ) ,C < 15 > 
C(10)=l. 
IBP  =  1 

X(5)=61.094 
PI=3. 14159? 
X(8)=IBP 
X(10)=0.C 
1    CALL    INTEG2(T,X,XD0T,C) 
X(10)=T-C. 097687 
XMAX=C(l)*SIN(PI/3.) 
THA=2.C*PI*60.C*T 
PHASEA=ABS<C(1)*SIN(THA)) 
PHASEB=ABS<C(l>*SIN<THA*2.0*PI/3.  )) 
PHASEC=ABS<C(l)*SIN(THA*PI/3.)) 
X(6)=AMAX1 (PHASEA.PHASEB,PHASEC) 
IF(PHASEA.GE.XMAXIlpH=l 
IF(PHASEB.GE.XMAXHPH=2 
IF(PHASEC.GE.XMAXUPH=3 

206  ERR0R=X(1)-C(5) 
X(8)=IPH 
VKl=C(6)*ERR0R 
IF(VK1.GE.C(7))VDZ=VK1-C(7) 
IF<  VKl.LT.Cm  )VDZ  =  0.0 

IF( VK1,LT.0.0)VDZ=VK1 

VK2=VDZ*C(8> 

GO   TO< IOC, 101, 102, 103, 104), IBP 

100  CONTINUE 
IF(VK2)2C7,207,208 

207  IBP=1 
VIN=X(6) 
X(5)=VIN 
X<7)=0. 

X<3)=(l./C(2)  )*(VIN-X(1>) 
XD0T(2)=X<3) 
IF<X(2).LT.-0.C1)  X(2)=-0.01 

xui*<i./c(3ii*x<2i-<i./<c(3)*cmn*xm 

XD0T(1)=X(4) 
VLAST=VIN 
X(8)=IBP 
GO  TO  1 

208  IF(VIN.LT.VLAST)GO  TO  207 
GO  TOCIOI.102,103) ,IPH 

101  VIN=PHASEA 
IF(VK2>207,207,105 

105  CONTINUE 
X(7)=l. 
X(5)=VIN 
IF(VIN-. 75)300,301,301 

300  VIN=0.0 
X(5)=VIN 
IBP  =  5 

GO  TO  302 

301  IBP=2 

302  X(3)=<1./C(2)I*<VIN-X(1)) 
IF(X<2).LT.-0.C1)X(3)=0.0 
XD0T(2)=X(3) 
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10? 
1C6 

30  3 


304 

3  J5 


103 
107 

306 


3C7 
30  3 


104 


309 
310 


IP(X( 
X(4)  = 

XDGT( 
VUST 
X(8)  = 

GO  TO 
VIN=P 
IF(  VK 

CONTI 

X<7)  = 

X<5)  = 
IF(VI 
VIN=0 
X(  5  >  = 

IBP  =  5 
GO  TO 
IRP  =  3 
X(3)  = 
IF(X( 
XDOT< 
IF(X( 
X(4)  = 
XDGT( 
VLAST 
X(8)  = 
GO  TO 
VIN  =  P 
IF(VK 
CONTI 
X(7)  = 
X(5)  = 
IFCVI 
VIN  =  C 
X(5)  = 
IBP  =  5 
GO  TO 
IBP  =  4 
X  (  3  >  = 
IF(X( 
XDOT( 
IF(X( 
X(4>  = 
XDOT( 
VLAST 
X<8)  = 
GO  TO 
VIN=C 
X(7)  = 
X(5)  = 
X<3)  = 
IF(X( 
XDOT( 
IF(X( 
X(4>  = 
XDOT( 
IBP=5 
IF(  VK 
IBP=1 
CONTI 
X(8)  = 
GO  TO 
END 


2).LT.-0.C1)  X(2)=-n.^l 
<1./C(3))*X<2)-(1./(C(3)*C(4)))*X(1) 

1)=X(M 

=  VIN 

IBP 
i 

HASEB 

2)207,207,106 

NUE 

2. 

VIN 

N-. 75)303, 304, 304 

•  0 

VIN 

3C5 

<1./C<2) IMVIN-X(l) ) 

2).LT.-0.C1)X(3»=0.0 

2)=X(3) 

?).LT.-G.fl)  X(2)=-0.01 

<1./C<3)MX<2>-(1./<C(3)»C(4)))*X(1) 

1)=X(4) 

=  VIN 
IBP 

1 
HASFC 

?) 207,207,107 
NUE 
3. 
VIN 

N-. 75) 306,307,307 
.0 
VIN 

308 

(1./CC2) )* (VIN-X(l) ) 

2>.LT.-o.mxm=o.G 

2)=X( 3) 

2).LT.-0.C1)     X(2)=-0.C1 

(l./C<  3))*X<2)-(1./(C(3)*C(4) )  )*X( 1) 

1)=X(4) 

=  VIN 

IBP 

1 
.0 
4. 
VIN 

U./CI2) )»(VIN-X(1)  ) 
2).LT.-0.(  1 IX<3)=0.0 
2)=X(3) 
2I.LT.-0.C1)     X<21=-0.01 

a,/c(3)i*x(2i-(i./(ci3i*cK)n*xm 

1)=X(4) 

2)309,309,310 

NUE 

IBP 

1 


81 


INPUT  DATA  RECORD 


ORDER  OF  EQUATIONS  =   2 

INITIAL  TIME  =  0.9769E-O1 

FINAL  TIME  =  0.2000E  OC 

STEP  SIZE  =  0.2222E-04 

THE  NON-ZERO  CONSTANTSt  CM),  ARE 

C(  1)  =  0.650CE  02 

C(  2)  =  0.5GCCE-G3 

C(  3)  =  0.10COE  CO 

C(  4)  =  0.320CE  CI 

C(  5)  =  0.600CE  02 

C<  6)  =  0.1C0CE  01 

C(  8)  =  0.1C00E  CI 

THE  NON-ZERO  INITIAL  CONDITIONS  ARE 

X(  1)  =  Q.6006E  02 

X(  21  =  0.3786E  02 


THE  COLUMN  HEADINGS  AND  THE  CORRESPONDING  VARIABLES  ARE 

TIME  X<  0) 

VIN  X(  5) 

VOUT  X(  1) 
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HYDE,W.H.   B16   P=3.20HMS,   VREF=60.0 VOLTS 


TIME  VIN  VOUT 


0.97687E- 

-01 

0.61094E 

02 

0.6H059E 

02 

C.9813GE- 

-01 

0.63935E 

02 

C.60148E 

02 

0.98574E- 

-01 

0.64994E 

02 

0.60251E 

02 

C.99C17E- 

-01 

0.6424CE 

02 

C.60372E 

02 

C.99461E- 

-01 

G.61695E 

C2 

0.60507E 

02 

L.999C4E- 

-01 

0.57430E 

02 

0.60647E 

02 

G.10O35E 

CO 

0.51563E 

02 

C.60773E 

02 

C.10C79E 

CO 

0.44258E 

02 

G.60862E 

02 

0.10123E 

00 

0.3572CE 

02 

0.60885E 

02 

C.10168E 

00 

0.26186E 

02 

C.60819E 

02 

0.10212E 

00 

0.15921E 

02 

C.60734E 

02 

O.10256E 

00 

0.52118E 

0! 

C.6065CE 

02 

0. 103C1E 

00 

0.0 

0.60566E 

02 

G.10345E 

00 

0.0 

0.60481E 

C2 

0.10390E 

00 

0.0 

C.60397F 

02 

C.10434E 

00 

o.c 

0.60313E 

02 

0.10478E 

00 

0.0 

0.60229E 

02 

C.1C523E 

00 

0.0 

0.60146E 

02 

C.10567E 

00 

0.0 

C.60062E 

02 

C.10611E 

00 

0.61829E 

02 

0.59978E 

02 

0.1C656E 

00 

J.643C4E 

02 

0.59901E 

02 

0.107C0E 

00 

0.64986E 

02 

0.59840E 

02 

C.1C74  4E 

00 

0.63856E 

02 

0.59799E 

02 

C.10789E 

00 

0.60945E 

0? 

C. S9773E 

02 

0.10833E 

00 

0.56335E 

02 

0.59751E 

02 

C.10877E 

00 

0.60885E 

02 

0.59723E 

02 

0.10922E 

00 

0.63823E 

02 

0.59698E 

02 

C.10966E 

CO 

0.649 82 E 

02 

0.59690E 

02 

C.11C10E 

00 

0.64329E 

C2 

0.59701E 

02 

C.11055E 

00 

0.61882E 

02 

0.59730E 

02 

0.11099E 

00 

0.57710E 

02 

0.59768E 

0  2 

0.11143E 

00 

0.59822E 

.02 

C.59799E 

02 

0.11188E 

00 

0.63218E 

02 

C.59829E 

02 

C.11232E 

CO 

0.64852E 

02 

C.59873E 

02 

0.11276E 

00 

0.64677E 

02 

C.59936E 

02 

0.  11321E 

00 

0.62699E 

02 

C.60016E 

02 

0.11365E 

00 

0.58973E 

02 

C.60107E 

02 

C.11409E 

00 

0.536C2E 

02 

0.60192E 

02 

C.11454E 

00 

0.46737E 

02 

0.60251E 

02 

0.11498E 

00 

0.38568E 

02 

0.60255E 

02 

G.11543E 

00 

0.29324E 

02 

0.60186E 

02 

0.11587E 

00 

0.19262E 

02 

C.60102E 

02 

0.11631E 

00 

0.86633E 

01 

0.60019E 

C2 

C.11676E 

00 

0.57348E 

02 

0.59935E 

C2 

0.11720E 

00 

0.61640E 

02 

0.59852E 

02 

C.11764E 

00 

0.64213E 

02 

0.59769E 

02 

0.11809E 

00 

0.64996E 

02 

C.59700E 

02 

G.11853E 

00 

0.63966E 

02 

0.59651E 

02 

0.11897E 

00 

0.61152E 

02 

0.59619E 

C2 

0.11942E 

00 

0.56634E 

02 

0.59592E 

02 
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HYDE,W.H.   B16   P=3.20HMS,   VREF=60.D VOLTS 


TIME 


VIN 


VOUT 


0.11986E 

00 

0.60671E 

02 

0.59559E 

02 

G.12C30E 

00 

0.637C6E 

02 

0.59529E 

02 

C.12075E 

GO 

0.64965E 

02 

C.59516E 

02 

0.12119E 

00 

0.64413E 

02 

C.59524E 

02 

G.12163E 

00 

0.62065E 

C2 

C.59550E 

02 

C.122C8E 

00 

0.57985E 

02 

0.59586E 

C2 

C.12252E 

00 

0.59583E 

02 

0.59617E 

02 

0.12296E 

00 

0.63C75E 

02 

C.59647E 

02 

0.12341E 

00 

0.64808E 

02 

0.59690E 

02 

0.12385E 

00 

G.64735E 

02 

C.59753E 

02 

C.12429E 

00 

0.62856E 

02 

0.59835E 

02 

0.12474E 

00 

0.59225E 

02 

0.59928E 

02 

0.12518E 

00 

0.58379E 

02 

0.60018E 

02 

0.12562E 

00 

0.62321E 

02 

C.60103E 

02 

0.126C7E 

00 

0.64525E 

02 

0.60195E 

02 

0. 12651E 

00 

0. 64930 E 

02 

C.60304E 

02 

0.12696E 

CO 

0.63525E 

02 

0.60430E 

02 

C.12740E 

00 

0.60348E 

C2 

0.60568E 

02 

0.  12784E 

00 

0.55489E 

02 

G.60704E 

02 

0.12829E 

00 

0.49C82E 

02 

C.60819E 

02 

0.12873E 

00 

0.41307E 

02 

0.60887E 

02 

C.12917E 

00 

0.32380E 

02 

C.60878E 

02 

C.12962E 

00 

0.22549E 

02 

0.6G796E 

02 

0.13006E 

00 

0.12C91E 

02 

0.60711E 

02 

0.13050E 

00 

0.12949E 

01 

0.60627E 

02 

0.13G95E 

00 

0.0 

0.60542E 

02 

0.13139E 

00 

0.0 

0.60458E 

02 

0.13183E 

00 

0*0 

0.60374E 

02 

0.13228E 

00 

o.c 

0.60290E 

02 

0.13272E 

00 

0.0 

0.60206E 

02 

C.13316E 

00 

0.0 

0.60122E 

02 

0.13361E 

00 

0.0 

0.60039E 

02 

0.13405E 

00 

0.62926E 

02 

0.59956E 

02 

0.13449E 

00 

0.64759E 

02 

0.59883E 

02 

0.13494E 

00 

0.64786E 

02 

0.59828E 

02 

0.13538E 

00 

0.63C07E 

02 

0.59793E 

02 

0.13582E 

00 

0.594  71E 

02 

0.5977CE 

02 

0.13627E 

00 

0.58111E 

02 

0.59746E 

02 

0.13671E 

00 

0.62147E 

02 

C.59716E 

02 

0.13715E 

00 

0.64450E 

02 

0.59696E 

02 

0.13760E 

00 

0.64956E 

02 

C.59693E 

02 

0.  138C4E 

00 

0.63650E 

02 

0.59711E 

02 

0.13849E 

00 

0.605 70E 

02 

0.59744E 

02 

0.13893E 

00 

0.56770E 

02 

0.59780E 

02 

0.1393  7E 

.88 

0.61246E 

02 

0.59808E 

02 

C.13982E 

0.64015E 

02 

0.59842E 

02 

0.14C26E 

00 

0.64998E 

02 

0.59891E 

02 

0.14070E 

00 

0.64170E 

02 

0.59960E 

02 

0.14115E 
0.14159E 

CO 

0.61552E 
0.57217E 

02 

0.60045E 

02 

00 

02 

0.60135E 

02 
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TIME 


VIN 


VOUT 


°« 

14203E 

00 

0.51287E 

02 

c  • 

14248E 

00 

0.43927E 

02 

0. 

14292E 

CO 

0.35343E 

02 

0. 

14336E 

00 

0.25772E 

02 

0. 

14381E 

C  ) 

0.15483E 

C2 

0. 

14425E 

CO 

J.58522E 

02 

c. 

14469E 

OC 

0.59090E 

L2 

0. 

14514E 

00 

0.62772E 

r? 

0, 

14558E 

00 

0 .64704  6 

02 

c. 

146C2F 

00 

C.64833E 

02 

0.14647E  00 

C.14691E  00 

0. 14735E  00 

L.14780E  CO 

C.14824E  00 

C.14868E  CO 

C. 14913E  00 

C.14957E  CO 

0.15OO2E  Oo 

0.15C46E  00 

G.15CQ0E  00 
C.15135E  00 
C.15179F  CO 
C.15223E  00 
0.15268E  CC 
C.15312E  00 
0.15356E  CO 
C.154C1E  00 
0.15445E  00 
C. 15489E  00 


0.15 
0.15 
0.15 
C.15 
0.  15 
0.15 
0.15 
0.15 
C.15 
0.15 


534E 
578E 
622E 
66  7E 
711E 
755E 
80  CE 
844E 
888E 
933E 


CO 
CO 
CO 
CO 
00 
00 
00 
00 
00 
CO 


0.15977E  CO 
C.16021E  CO 
0.16066E  00 
C.16110E  CO 
0.16154E  00 
C.16199E  00 
0.16243E  00 
0.16288E  CO 
0.16332E  00 
0.16376E  00 


0.63153E  02 
0.59712E  02 
0.57837E  02 
0.61967E  02 
0.64368E  02 
0.64975E  C2 
0.63770E  02 
0.60787E  C2 
0.56473E  02 
0.61C41E    0? 


C. 63907 
0.64991 
1.64263 
0.61743 
0.5750J 
0.59996 
0.63321 
0.64881 
0.64631 
0.62579 


02 
02 
02 
02 
02 
02 
02 
C2 
02 
02 


0.5  8782E  02 
0.53346E  02 
0.46423E  02 
0.38205E  02 
0.28921E  02 
0.18831E  02 
0.82169E  01 
0.0 
0.0 
0.0 

O.C 

0.0 

0.0 

0.0 

0.0 

0.63794E  C2 

0.64978E  02 

0.64351E  02 

0.61929E  C2 

0.57781E  02 


C. 60213 
0.60255 
C. 60232 
0.60150 
C.6CC66 
0.59983 
0.598^9 
0.59836 
0.59736 
0.59674 


0.59656 
0.597r6 
C. 59777 
0. 5^864 
0.59958 
0.60045 
C.6C132 
0.6023C 
0. 60346 
C. 60477 


0.606 
0.607 
0.6C8 
C.608 
C.608 
0.607 
C.6G6 
C.605 
C.605 
0.6C4 


0.60344 
0.6C2  60 
C.6C176 
0.60092 
0.60008 
0.59927 
0.59859 
0.59811 
0.59780 
0.59757 


02 
02 

02 
02 
G2 
02 
02 
02 
0  2 
02 


0.596321 

0.596031 

C.59575E    02 

0.59541 

0.5<?5]7| 

C.59511 

0.595261 

C. 595571 

0.595931 

C.59622I 


0  2 
02 


G2 
02 
02 
02 
02 
02 
02 


02 
02 
02 
02 
02 
02 
02 
02 
02 
02 


17E  C2 
48E  02 
50E  02 
94E  02 
5CE  02 
65E  02 
81E  02 
96E  02 
12E  02 
28E  02 


02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
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HYDE,W.H.   B16   R=3.20HMS,   VREF=60.0 VOLTS 


TIME 


VIN 


VOUT 


0.16421E  00 
0.16465E  00 
0.165G9E  00 
0.16554E  00 
0.16598E  00 
0.16642E  00 
0.16687E  00 
0.16731E  00 
0.16775E  00 
0.16820E  00 

0.16864E  00 
©•16908E  00 
0.16953E  00 
0.16997E  00 
0.17G41E  00 
C.17086E  00 
C.17130E  CO 
0.17174E  CO 
0.17219E  00 
C.17263E  00 

0.17307E  00 
C.17352E  00 
C.17396E  00 
0.17441E  00 
G.17485E  00 
0.17529E  00 
0.17574E  00 
0,17618E  00 
0.17662E  00 
C.17707E  00 

C.17751E  00 
C.17795E  00 
0.17840E  00 
0.17884E  00 
C.17928E  00 
0.17973E  00 
0.18017E  00 
0.18C61E  00 
C.18106E  00 
C.18150E  00 

0.18194E  00 
0.18239E  00 
0.18283E  00 
0.18327E  00 
0, 18372E  00 
0.18416E  00 
0.18460E  00 
0.18505E  00 
0.18549E  00 
0.18594E  00 


0.59761E  02 
0.63182E  02 
0.64841E  C2 
0.64693E  C2 
0.62740E  02 
0.59037E  02 
0.58576E  02 
0.62448E  C2 
0.64578E  02 
0.649C8E  02 

0.63428E  02 
0.60179E  02 
0.55253E  02 
0.48786E  02 
0.40957E  02 
0.31988E  02 
0.22126E  02 
0.11648E  02 
0.567C9E  C2 
0.6C615E  02 

0.63675E  02 
0.64960E  02 
0.64433E  C2 
0.62110E  02 
0.53055E  02 
0.59521E  02 
0.63038E  02 
0.64796E  02 
0.64748E  02 
0.62895E  02 

0.59288E  02 
0.58312E  02 
0.62277E  02 
0.64507E  02 
0.64937E  02 
0.63557E  02 
0.60405E  02 
0.56988E  02 
0.61396E  02 
0.64091E  02 

0.65000E  02 
0.64096E  02 
0.61405E  02 
0.57002E  02 
0.51009E  02 
0.43594E  02 
0.34963E  02 
0.25358E  02 
0.15045E  02 
0.43128E  01 


C.59729E  G2 
0.59700E  02 
0.59685E  02 
0.59689E  02 
C.59713E  02 
0.59748E  02 
0.5978CE  02 
0.59807E  02 
0.59845E  02 
0.59901E  02 

0.59976E  02 
0.6C064E  02 
C.6015IE  02 
G.60219E  02 
0.60241E  02 
0.6C188E  02 
0.6G1C4E  02 
0.60021E  02 
0.59937E  02 
0.59854E  02 

0.59771E  02 
0.59695E  02 
0.59639E  02 
0.59602E  02 
0.59574E  02 
0.59541E  02 
0.59509E  02 
0.59489E  02 
0.59490E  02 
0.59512E  02 

0.59546E  02 
0.59578E  02 
0.596C7E  02 
0.59646E  02 
0.59703E  02 
0.59780E  02 
0.59872E  02 
0.59965E  02 
0.60050E  02 
C.60140E  02 

0.60244E  02 
0.60367E  02 
0.60504E  02 
0.60643E  02 
0.60768E  02 
0.60853E  02 
0,60869E  02 
0.60799E  02 
0.60715E  02 
0.60630E  02 


86 


HYOE,W.H.   B16   P=3.20HMS,   VREF=60.D VOLTS 


TIME 


VIN 


VOUT 


0.18638E 

00 

0.0 

0.60546E 

02 

0.18682E 

00 

3.0 

0.60462E 

02 

0.18727E 

00 

3.0 

C.60377E 

C2 

0.18771E 

00 

0.0 

0.6C294E 

02 

0.18815E 

00 

0.0 

C.6G210E 

C2 

C.18860E 

00 

0.0 

0.60126E 

02 

C. 189C4E 

00 

0.0 

C.60042E 

02 

0.18948E 

00 

0.62101E 

02 

C.59959E 

02 

0.18993E 

CO 

0.64430E 

02 

0.59881E 

02 

G.19C37E 

CO 

0.64961E 

C2 

C.59821E 

02 

0.19081E 

00 

0.63681E 

02 

0.59781E 

02 

C.19126E 

00 

0.60626E 

C2 

0.59755E 

02 

0.19170E 

00 

0.56695E 

02 

0.59732E 

02 

C.19214E 

00 

0.61195E 

02 

0.597C3E 

02 

C.19259E 

CO 

0.63988E 

02 

0.59679E 

C2 

C.19303E 

00 

0.64997E 

C2 

0.59671E 

02 

0.19347E 

00 

0.64194E 

02 

0.59684E 

02 

0.19392E 

00 

0.61601E 

02 

C.59714E 

02 

0.19436E 

CO 

0.57290E 

02 

0.59751E 

02 

0.19480E 

00 

0.60169E 

02 

C.59781E 

02 

C.19525E 

00 

0.63422E 

02 

0.59813E 

02 

G.19569E 

00 

0.64907E 

02 

0.59858E 

02 

C.19613E 

00 

0.64581E 

02 

0.59922E 

02 

0.19658E 

00 

J.62456E 

02 

0.60004E 

02 

0.19702E 

00 

0.58588E 

C2 

0.60095E 

C2 
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INPUT  DATA  RECORD 


ORDER  OF  EQUATIONS  =   2 
INITIAL  TIME        *  0.9769E-01 
FINAL  TIME  =  0.2000E  00 

STEP  SIZE  =  0.2222E-04 

THE 


NON-ZERO 

CONSTANTS,  C(I),  ARE 
C(  1)  =  0.650CE  02 
C(  2)  =  0.5000E-03 
C<  3)  =  O.IOOOE  00 
C(  4)  =  0.1600E  01 
C(  5)  =  0.6000E  02 
C(  6)  =  O.IOOOE  01 
C<  8)  =  O.IOOOE  01 

NON-ZERO 

INITIAL  CONDITIONS  ARE 
X(  1)  =  0.5897E  02 
X(  2)  =  0.4501E  02 

COLUMN  HEADINGS  AND  THE  CORRESPONDING  VARIABLES 

TIME         X(  01 
VIN          X(  5) 
VOUT         X<  1> 

THE 


THE  COLUMN  HEADINGS  AND  THE  CORRESPONDING  VARIABLES  ARE 
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HYDE,W.H.   Bl  R=1.6  OHMS,   VREF=  63.0  VOLTS. 


TIME 


VIN 


VOUT 


0.97687E- 

-01 

0.61094E 

02 

0.58968F 

02 

C.98130E- 

-CI 

0.63935E 

02 

0.59010E 

02 

0.98574E- 

-01 

0.649 94E 

02 

0.  5  9072E 

02 

0.99017E- 

-01 

0.6424CE 

02 

0.59155E 

02 

C.99461E- 

-01 

0.61695E 

02 

0.5925RE 

02 

0.99904E- 

-01 

0.57430E 

02 

0.5  93  70E 

02 

C.10035E 

00 

0.60055E 

02 

0.59476E 

02 

0.1C079E 

00 

0.63356E 

02 

0.59584E 

02 

0.10123E 

00 

0.64890E 

02 

C.59706E 

C2 

C.10168E 

00 

0.64614E 

02 

0.59848E 

02 

0.10212E 

00 

0.62537E 

02 

0.600G8E 

02 

0.10256E 

00 

0.58716E 

02 

0.60176E 

02 

0.103C1E 

00 

0.53258E 

02 

0.60338E 

02 

0.10345E 

00 

0.46315E 

02 

0.60470E 

02 

0.10390E 

CO 

0.38080E 

02 

0.60546E 

02 

C.10434E 

00 

0.28784E 

02 

0.60533E 

02 

0.10478E 

00 

0.18685E 

02 

0.60402E 

02 

0.1O523E 

00 

0.80639E 

01 

0.60234E 

02 

0.1056  7E 

00 

0.0 

0.60067E 

02 

C.10611E 

00 

0.61829E 

02 

0.59900E 

02 

0.10656E 

00 

0.64304E 

02 

0.59739E 

02 

C.107C0E 

00 

0.64986E 

02 

0.59596E 

02 

0.10744E 

00 

0.63856E 

02 

0.59474E 

02 

0.10789E 

00 

0.60945E 

02 

0.59369E 

02 

0.10833E 

00 

0.56335E 

02 

0.59270E 

02 

0.1C877E 

00 

0.60885E 

02 

0.59166E 

02 

O.10922E 

00 

0.63823E 

02 

0.59069E 

02 

0.10966E 

00 

0.64982E 

02 

C.58990E 

02 

O.llOlOE 

00 

0.64329E 

02 

0.58934E 

02 

0.11055E 

00 

0.61882E 

02 

0.58899E 

02 

0.11099E 

00 

0.57710E 

02 

C.58876E 

02 

0.11143E 

00 

0.59822E 

02 

0.58850E 

02 

0.11188E 

00 

0.63218E 

02 

0.58827E 

02 

0.11232E 

00 

0.64852E 

02 

0.58822E 

02 

0.11276E 

00 

0.64677E 

02 

0.58840E 

02 

C.U321E 

00 

0.62699E 

02 

0.58880E 

02 

0.11365E 

00 

0.58973E 

02 

0.58934E 

02 

0.11409E 

00 

0.58642E 

02 

0.58989E 

02 

0.11454E 

00 

0.62490E 

02 

0.59042E 

02 

0.11498E 

00 

0.64595E 

02 

0.59109E 

02 

0.U543E 

00 

0.64900E 

02 

0.59196E 

02 

C.11587E 

00 

0.63394E 

02 

0.59305E 

02 

0.11631E 

00 

0.60121E 

02 

0.59429E 

02 

0.11676E 

00 

0.57348E 

02 

0.59555E 

02 

0.11720E 

00 

0.61640E 

02 

0.59675E 

02 

0.11764E 

00 

0.64213E 

02 

0.59802E 

02 

0.118C9E 

00 

0.64996E 

02 

C.59946E 

02 

0.11853E 

00 

0.63966E 

02 

0.60108E 

02 

0.11897E 

00 

0.61152E 

02 

0.60285E 

02 

0.11942E 

00 

0.56634E 

02 

0.60464E 

02 
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HYDE,W.H.   Bl  R=1.6  OHMS,   VREF=  60.0  VOLTS. 


TIME 


VIN 


VOUT 


0.11986E  00 
0.12030E  00 
0.12075E  00 
0.12119E  00 
0.12163E  00 
0.12208E  00 
0.12252E  00 
0.12296E  00 
0.12341E  00 
0.12385E  00 

C.12429E  00 
0.12474E  00 
0.12518E  00 
0.12562E  00 
0.12607E  00 
0.12651E  00 
0.12696E  00 
0.12740E  00 
0.12784E  00 
0.12829E  00 

0.12873E  00 
0.12917E  00 
0.12962E  00 
0.13006E  00 
0.13C50E  00 
0.13095E  00 
0.13139E  00 
0.13183E  00 
0.13228E  00 
0.13272E  00 

0.13316E  00 
0.13361E  00 
C.13405E  00 
0.13449E  00 
0.13494E  00 
0.13538E  00 
0.13582E  00 
0.13627E  00 
0.13671E  00 
0.13715E  00 

0.13760E  00 
0.13  804E  00 
0.13849E  00 
0.13893E  00 
0.13937E  00 
0.13982E  00 
0.14026E  00 
0.14070E  00 
0.14115E  00 
0.14159E  00 


0.50536E  02 
0.43029E  02 
0.34322E  02 
0.24658E  02 
0.14307E  02 
0.35561E  01 
0.0 
0.0 

0.64808E  02 
0.64735E  02 

0.62856E  02 
0.59225E  02 
0.58379E  02 
0.62321E  02 
0.64525E  02 
0.64930E  02 
0.63525E  02 
0.60348E  02 
0.57061E  02 
0.61446E  02 

0.64117E  02 
0.65000E  02 
0.64070E  02 
0.61354E  02 
0.56928E  02 
0.60451E  02 
0.63583E  02 
0.64943E  02 
0.64491E  02 
0.62241E    02 

0.58256E  02 
0.59339E  02 
0.62926E  02 
0.64759E  02 
0.64786E  02 
0.63007E  02 
0.59471E  02 
0.54277E  02 
0.47569E  02 
0.39534E    02 

0.30397E  02 
0.20413E  02 
0.98602E  01 
0.0 

0.61246E  02 
0.64015E  02 
0.64998E  02 
0.64170E  02 
0.61552E  02 
0.57217E   02 


0.60627E  02 
0.60749E  02 
0.60800E  02 
0.60747E  02 
0.605  87E  02 
0.60419E  02 
0.60251E  02 
0.60083E  02 
0.59919E  02 
0.59771E  02 

0.59642E  02 
0.59526E  02 
0.59409E  02 
0.59289E  02 
0.59181E  02 
0.59093E  02 
0.59029E  02 
0.58981E  02 
0.58939E  02 
0.58893E  02 

0.58857E  02 
0.58841E  02 
0.58849E  02 
0.58877E  02 
0.5891AE  02 
0.58947E  02 
0.58986E  02 
0.59042E  02 
0.59121E  02 
0.59220E  02 

0.59330E  02 
0. 59436E  02 
0.59542E  02 
C.59661E  02 
0.59798E  02 
0.59954E  02 
0.60122E  02 
0.60285E  02 
0.60424E  02 
0.60512E  02 

0.60516E  02 
0.60402E  02 
0.60234E  02 
0.60067E  02 
0.59900E  02 
0.59737E  02 
0.59590E  02 
0.59465E  02 
0.59358E  02 
0.59259E  02 
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HYDEfW.H.   Bl  R  =  1.6  OHMS,   VREF=  60,0  VOLTS. 


TIME 


VIN 


VOUT 


0.14203E 

00 

0.60226E 

02 

0.59156E 

02 

0.14248E 

00 

0.63455E 

02 

0.59057E 

02 

0.14292E 

00 

0.64915E 

02 

0. 58975E 

02 

0.14336E 

00 

0.64564E 

02 

0.58916E 

02 

0.14381E 

CO 

0.62413E 

02 

0.58878E 

02 

0.14425E 

CO 

0.58522E 

02 

0.58855E 

C2 

0.14469E 

00 

0.59090E 

02 

0.58830E 

02 

0.14514E 

00 

0.62772E 

02 

0.58807E 

02 

0.14558E 

00 

0.64704E 

02 

0.58799E 

02 

0. 14602E 

CO 

0.64833E 

02 

0.58814E 

02 

0.14647E 

00 

0.63153E 

02 

0.58852E 

02 

0.14691E 

00 

0.59712E 

02 

0.58906E 

02 

0.14735E 

00 

0.57837E 

02 

0.58963E 

02 

0.14780E 

00 

0.61967E 

02 

0.59017E 

02 

0.14824E 

00 

0.64368E 

02 

C.59082E 

02 

C.14868E 

00 

0.64975E 

02 

0.59167E 

02 

0.14913E 

00 

0.63770E 

02 

0.59274E 

02 

0.14957E 

00 

0.60787E 

02 

0.59398E 

02 

0.15002E 

00 

0,564 73E 

02 

C.59527E 

C2 

0.15046E 

00 

0.61041E 

02 

0.59649E 

02 

0.15090E 

00 

0.63907E 

02 

0.59776E 

02 

0.15135E 

00 

0.64991E 

02 

0.59918E 

02 

0.15179E 

00 

0.64263E 

02 

0.60079E 

02 

0.15223E 

CO 

0.61743E 

C2 

0.60255E 

02 

0.15268E 

00 

0.57501E 

02 

0.60437E 

02 

0.15312E 

00 

0.51656E 

02 

0.60605E 

02 

0.15356E 

00 

0.44371E 

02 

0.60738E 

02 

0.15401E 

00 

0.35848E 

02 

0.60805E 

02 

0.15445E 

no 

0.26326E 

02 

C.60773E 

02 

0.15489E 

00 

0.  160  70E 

02 

0.60625E 

02 

0.15534E 

00 

0.53645E 

01 

0.60456E 

02 

0.15578E 
0.15622E 

00 

0.0 

0.60288E 

02 

00 

0.0 

0.60121E 

02 

0.15667E 

00 

0.64644E 

02 

0.59954E 

02 

0.15711E 

00 

0.64873E 

02 

0.59801E 

02 

0.15755E 

CO 

0.63293E 

02 

0.59668E 

02 

0.15800E 

00 

0.5994QE 

02 

0.59549E 

02 

0.15844E 

00 

0.57559E 

02 

0.59431E 

02 

0.15888E 

00 

0.61782E 

02 

0.59309E 

02 

0.15933E 

00 

0.64282E 

02 

0.59197E 

02 

0.15977E 

00 

0.64989E 

02 

C.59104E 

02 

0.16021E 

00 

0.63884E 

02 

0.59034E 

02 

O.16066E 

00 

0.60998E 

02 

0.58983E 

02 

0.16110E 

00 

0.56411E 

02 

0.58939E 

02 

0.16154E 

00 

0.60831E 

02 

0.58891E 

02 

0.16199E 

00 

0.63794E 

02 

0.58851E 

02 

0.16243E 

00 

0.64978E 

Q2 

0.58830E 

02 

0.16288E 
0.I6332E 

00 

0.64351E 

02 

0.58833E 

02 

CO 

0.61929E 
0.57781E 

02 

0.58856E 

02 

0.16376E 

00 

02 

0.58892E 

02 
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HYDE,W.H.   Bl  R=1.6  OHMSt   VREF*  60.0  VOLTS 


TIME 


VIN 


VOUT 


0.16A21E  00 
0.16465E  00 
0.16509E  00 
0.16554E  00 
0.16598E  00 
0. 16642E  00 
0.16687E  00 
0.16731E  00 
0.16775E  00 
0.16820E  00 

0.16864E  00 
0.16908E  00 
0.16953E  00 
0.16997E  00 
0.17041E  00 
0.17C86E  00 
0.I7130E  00 
0.17174E  00 
0.17219E  00 
0.17263E  00 

0.17307E  00 
0.17352E  00 
0.17396E  00 
0.17441E  00 
0.17485E  00 
0.17529E  00 
0.17574E  00 
0.17618E  00 
0.17662E  00 
0.17707E  00 

0.17751E  00 
0.17795E  00 
0.17  840E  00 
0.17884E  00 
0.17928E  00 
0.L7973E  00 
0.18017E  00 
0.18061E  00 
0.18106E  00 
0.18150E  00 

0.18194E  00 
0.18239E  00 
0.18283E  00 
0.18327E  00 
0.18372E  00 
0.1 841 6E  00 
0.18460E  00 
0.18505E  00 
0.18549E  00 
0.18594E  00 


0.59761E  02 
0.63182E  02 
0.64841E  02 
0.64693E  02 
0.62740E  02 
0.59037E  02 
0.58576E  02 
0.62448E  02 
0.64578E  02 
0.64908E  02 

0.63428E  02 
0.60179E  02 
0.55253E  02 
0.48786E  02 
0.40957E  02 
0.31988E  02 
0.22126E  02 
0.U648E  02 
0.84376E  00 
0.60615E    02 

0.63675E  02 
0.64960E  02 
0.64433E  02 
0.62110E  02 
0.58055E  02 
0.59521E  02 
0.63038E  02 
0.64796E  02 
0.64748E  02 
0.62895E  02 

0.59288E  02 
0.58312E  02 
0.62277E  02 
0.64507E  02 
0.64937E  02 
0.63557E  02 
0.60405E  02 
0.56988E  02 
0.61396E  02 
0.64091E  02 

0.65000E  02 
0.64096E  02 
0.61405E  02 
0.57002E  02 
0.60394E  02 
0.63551E  02 
0.64936E  02 
0.64510E  02 
0.62285E  02 
0.58324E  02 


0.58924E  02 

0.58960E  02 

0.59012E  02 

0.59085E  02 

0.59181.E  02 

C.59289E  02 

0.59396E  02 

C.59500E  02 

0.59615E  02 

0.59749E  02 

0.59902E  02 

0.60068E  02 

0.6C233E  02 

0.60378E  02 

0.60476E  02 

0.60497E  02 

0.60405E  02 

0.60237E  02 

0.60070E  02 

0.59903E  02 

0. 59738E  02 

0.59586E  02 

0.59455E  02 

0.59344E  02 

0.59243E  02 

0.59140E  02 

0. 59038E  02 

0.58951E  02 

0.58887E  02 

0.58846E  02 

0.58821E  02 

0.58797E  02 

0.58772E  02 

0.58761E  02 

0.58772E  02 

0.58806E  02 

0.58859E  02 

0.58916E  02 

0.58971E  02 

•0.59034E  02 

0.59116E  02 

0.59221E  02 

0.59344E  02 

0.59474E  02 

0.59599E  02 

0.59725E  02 

0.59866E  02 

0.60026E  02 

0.60202E  02 

0.60385E  02 
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HYDE,W.H.   Bl  R=1.6  OHMS,   VREF  =  60,0  VOLTS 


TIME 


VIN 


VOUT 


0.18638E  00 
C.18682E  00 
0.18727E  00 
0.18771E  00 
0.18815E  00 
C.18860E  00 
0.18904E  00 
0.18948E  00 
C.18993E  00 
0.1903  7E  00 

0.19081E  00 
0.19126E  00 
C.19170E  00 
0.19214E  00 
0.19259E  00 
0.19303E  00 
0.19347E  00 
0.19392E  00 
0.19436E  00 
C.19480E  00 

0.19525E  00 
0.19569E  00 
0.19613E  00 
0.19658E  00 
0.19702E  00 


0.52736E  02 
0.45678E  02 
0.37346E  02 
0.27973E  02 
0.17820E  02 
0.71687E  01 
0.0 
0.0 

0.64430E  02 
0.64961E  02 

0.63681E  02 
0.60626E  02 
0.56695E  02 
0.61195E  02 
0.63988E  02 
0.64997E  02 
0.64194E  02 
0.61601E  02 
0.57290E  02 
0.60169E  02 

0.63422E  02 
0.64907E  02 
0.64581E  02 
0.62456E  02 
0.58588E  02 


0.6C560E  02 
0.60702E  C2 
0.60785E  02 
0.60774E  02 
0.60641E  02 
0.60473E  02 
0.60305E  02 
0.60138E  .02 
0.59971E  02 
0.59814E  02 

0.59677E  02 
0.59556E  02 
C.59439E  02 
0.59316E  02 
0.59201E  02 
0.59104E  02 
0.59030E  02 
0.58976E  02 
0.58932E  02 
0.58885E    02 

0.58843E  02 
0.58818E  02 
0.58817E  02 
0.58838E  02 
0.58872E  02 
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INPUT  DATA  RECORD 


ORDER  OF  EQUATIONS  =   2 
INITIAL  TIME        =  0.9769E-CI 
FINAL  TIME  =  0.2O0OE  00 

STEP  SIZE  =  0.2222E-04 

THE  NON-ZERO  CONSTANTSf  CCI),  ARE 

C(  1)  *  0.6500E  02 

C(  2)  =  0.5000E-03 

C(  31  =  0.1000E  00 

C(  4)  =  0.8000E  00 

C(  5)  =  0.6000E  02 

C(  6)  =  0.1000E  01 

C(  8)  =  0.1000E  01 

THE  NON-ZERO  INITIAL  CONDITIONS  ARE 
X(  1)  =  0.5709E  02 
X(  2)  =  0.1069E  03 


THE  COLUMN  HEADINGS  AND  THE  CORRESPONDING  VARIABLES  ARE 

TIME         X(  0) 
VIN  X(  5) 

VOUT         X(  1) 
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HYDE,W.H.   83   R  =  0.8  OHMS,  VREF=60.0 


TIME 


VIN 


VOUT 


0.97687E- 

•01 

0.61094E 

02 

0.57091E 

02 

C.9813CE- 

■01 

0.63935E 

02 

0.57258E 

02 

0.98574E- 

-01 

0.64994E 

02 

C.57450E 

02 

0.99017E- 

■01 

0.64240E 

02 

0.57670E 

02 

0.99461E- 

■01 

0.61695E 

02 

0.57915E 

02 

0.99904E- 

■01 

0.57430E 

02 

0.58172E 

02 

O.10C35E 

CO 

0.60055E 

02 

0.58427E 

C2 

0.10079E 

00 

0.63356E 

02 

C.58687E 

02 

0.10123E 

00 

0.64890E 

02 

0.58964E 

02 

0.10168E 

00 

0.64614E 

02 

0.59262E 

02 

0.10212E 

00 

0.62537E 

02 

0.59578E 

02 

0.10256E 

00 

0.58716E 

02 

0.59904E 

02 

0.10301E 

00 

0.53258E 

02 

0.60223E 

02 

0.10345E 

00 

0.46315E 

02 

0.60512E 

C2 

0.1039CE 

00 

0.38080E 

02 

0.60744E 

02 

0.10434E 

00 

0.28784E 

02 

0.60884E 

02 

0.10478E 

00 

0.18685E 

02 

0.60897E 

02 

0.10523E 

00 

0.80639E 

01 

0.60743E 

02 

Q.1C56  7E 

00 

0,0 

0.60424E 

02 

0.10611E 

00 

0.0 

0.60089E 

02 

0.10656E 

00 

0.64304E 

02 

0.59757E 

02 

0.10700E 

00 

0.64986E 

02 

0.59442E 

02 

0.10744E 

00 

0.63856E 

02 

0.59149E 

02 

0.10789E 

00 

0.60945E 

02 

0.58877E 

02 

0.10833E 

00 

0.56335E 

02 

0.58613E 

02 

0.10877E 

00 

0.60885E 

02 

0.58348E 

02 

0.1C922E 

00 

0.63823E 

02 

0. 58094E 

02 

0.10966E 

00 

0.64982E 

02 

0.57863E 

02 

0.11C10E 

00 

0.64329E 

02 

0.57661E 

02 

0.11055E 

00 

0.61882E 

02 

0.57486E 

02 

0.11099E 

00 

0.57710E 

02 

0.57329E 

02 

0.11143E 

00 

0.59822E 

02 

0.57176E 

02 

0.11188E 

oc 

0.63218E 

02 

0.57034E 

02 

0.11232E 

00 

0.64852E 

02 

0.56916E 

02 

0.11276E 

00 

0.646 77E 

02 

0. 56830E 

02 

0.11321E 

00 

0.62699E 

02 

0.56775E 

02 

0.11365E 

00 

0.58973E 

02 

0.56743E 

02 

0.11409E 

00 

0.58642E 

02 

0.56720E 

02 

0. 11454E 

00 

0.62490E 

02 

0.56705E 

02 

0.11498E 

00 

0.64595E 

02 

0.56712E 

02 

0.11543E 

00 

0.64900E 

02 

0. 56750E 

02 

0.11587E 

00 

0.63394E 

02 

0.56819E 

02 

0.11631E 

00 

0.60121E 

02 

0.56913E 

02 

0.11676E 

00 

0.57348E 

02 

0.57019E 

02 

0.11720E 

00 

0.61640E 

02 

0.57128E 

02 

0.11764E 

00 

0.64213E 

02 

0.57254E 

02 

0.11809E 

00 

0.64996E 

02 

0.57406E 

02 

0.11853E 

00 

0.63966E 

02 

0.57587E 

02 

0.11897E 

00 

0.61152E 

02 

0.57791E 

C2 

0.11942E 

00 

0.56634E 

02 

O.S8007E 

02 
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HYDE,W.H.   B3   R=0. 8  OHMSf  VREF=60.0 


TIME 


VIN 


VOUT 


0.11986E 

00 

0.60671E 

02 

0.58221E 

02 

0.12030E 

00 

0.63706E 

02 

C.58444E 

02 

0.12075E 

00 

0.64965E 

02 

0.58685E 

02 

0.12119E 

00 

0.64413E 

02 

0.58949E 

02 

C.12163E 

00 

0.62065E 

02 

C.59233E 

02 

0.12208E 

00 

0.57985E 

02 

0.59525E 

02 

0.12252E 

00 

0.59583E 

02 

0.59812E 

02 

0.12296E 

00 

0.63075E 

02 

0.60095E 

02 

0.12341E 

00 

0.64808E 

02 

0.60388E 

02 

0.12385E 

00 

0.64735E 

02 

0.60697E 

02 

0.12429E 

00 

0.62856E 

02 

0.61019E 

02 

0.12474E 

00 

0.59225E 

02 

0.61346E 

02 

0.12518E 

00 

0.53942E 

02 

0.61662E 

02 

0.12562E 

00 

0.47155E 

02 

0.61946E 

02 

0.12607E 

00 

0.39053E 

02 

0.62169E 

02 

0.12651E 

00 

0.29862E 

C2 

0.62299E 

02 

0.12696E 

00 

0.19839E 

02 

0.623C1E 

02 

0.1274CE 

00 

0.92624E 

01 

0.62136E 

02 

0.12784E 

00 

0.0 

0.61807E 

02 

0.12829E 

00 

0.0 

0.61465E 

02 

0.12873E 

00 

0.0 

0.61124E 

02 

0.12917E 

00 

0.0 

0.60785E 

02 

0.12962E 

00 

0.0 

0.60449E 

02 

0.13006E 

00 

0.0 

0.60113E 

02 

0.13050E 

00 

0.56928E 

02 

0.59780E 

02 

0.13C95E 

00 

0.60451E 

02 

0.59449E 

02 

0.13139E 

00 

0.63583E 

02 

0.59119E 

02 

0.13183E 

00 

0.64943E 

02 

0.58802E 

02 

0.13228E 

00 

0.64491E 

02 

0.58511E 

02 

0.13272E 

00 

0.62241E 

02 

0.58243E 

02 

0.13316E 

00 

0.58256E 

02 

0.57993E 

02 

0.13361E 

00 

0.59339E 

02 

C.57746E 

02 

0.13405E 

00 

0.62926E 

02 

0.57506E 

02 

0.13449E 

00 

0.64759E 

02 

0.57289E 

02 

0.13494E 

00 

0.64786E 

02 

0.57102E 

02 

0.13538E 

00 

0.63007E 

02 

0.56945E 

02 

0.13582E 

00 

0.59471E 

02 

0.56812E 

02 

0.13627E 

00 

0.58111E 

02 

0.56691E 

02 

0.13671E 

00 

0.62147E 

02 

0.56577E 

02 

0.13715E 

00 

0.64450E 

02 

0.56485E 

02 

0.13760E 
0.13804E 

00 

0.64956E 

02 

0.56424E 

02 

00 

0.63650E 

02 

0.56397E 

02 

0.13849E 

00 

0.60570E 

02 

0.56398E 

02 

0.13893E 

00 

0.56770E 

02 

0. 56415E 

02 

0,1393 7E 

00 

0.61246E 

02 

0.56438E 

02 

0.13982E 

00 

0.64015E 

02 

0.56479E 

02 

0.14026E 

00 

0.64998E 

02 

0.56549E 

02 

0.14070E 

00 

0.64170E 

02 

0.56651E 

02 

0.14115E 

00 

0.61552E 

02 

0.56782E 

02 

0.14159E 

00 

0.57217E 

02 

0.56930E 

02 
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HYDE.W.H.   B3   R=0.8  OHMS,  VREF=60.0 


TIME 


VIN 


VOUT 


0.14203E  00 
0.14248E  00 
0.14292E  00 
0.14336E  00 
0.14381E  00 
0.14425E  00 
0.14469E  00 
0.14514E  00 
0.14558E  00 
0.14602E  00 

0.14647E  00 
0.14691E  00 
0.14735E  00 
0.14780E  00 
0.14824E  00 
0.14868E  00 
0.14913E  00 
0.14957E  00 
G.15002E  00 
0.15046E  00 

0.15090E  00 
0.15135E  00 
0.15179E  00 
0.15223E  00 
0.15268E  00 
0.15312E  00 
0.15356E  00 
0.15401E  00 
0.15445E  00 
0.15489E  00 

0.15534E  00 
0.15578E  00 
0.15622E  00 
0.15667E  00 
0.15711E  00 
0.15755E  00 
0.15800E  00 
C.15844E  00 
0.15888E  00 
0.15933E  00 

0.15977E  00 
0.16021E  00 
0.16066E  00 
0.16110E  00 
0.16154E  00 
0.16199E  00 
0.16243E  00 
0.16288E  00 
0.16332E  00 
0.16376E    00 


0.60226E 

02 

0.57082E 

02 

0.63455E 

02 

0. 57245E 

02 

0.64915E 

02 

0.57431E 

02 

0.64564E 

02 

C.57645E 

02 

0.62413E 

02 

C.57884E 

02 

0.58522E 

02 

0.58140E 

02 

0.59090E 

02 

0.58396E 

02 

0.62772E 

02 

0.58653E 

02 

0.64704E 

02 

0.58925E 

02 

0.64833E 

02 

0.59217E 

02 

0.63153E 

02 

0.59529E 

02 

0.59712E 

02 

0.59853E 

02 

0.54607E 

02 

0.60175E 

02 

0.47978E 

02 

0.60472E 

02 

0.40C12E 

02 

0.60717E 

02 

0.30930E 

02 

0.6088CE 

02 

0.20986E 

02 

0.60923E 

02 

0.10456E 

02 

0*60809E 

02 

0.0 

0.60513E 

02 

0.0 

0.60178E 

02 

0.63907E 

02 

0.59845E 

02 

0.64991E 

02 

0.59528E 

02 

0.64263E 

02 

0.59233E 

02 

0.61743E 

02 

0.58960E 

02 

0.57501E 

02 

0.58698E 

02 

0.59996E 

02 

0.58436E 

02 

0.63321E 

02 

0.58180E 

02 

0.64881E 

02 

0.57946E 

02 

0.64631E 

02 

0.57740E 

02 

0.62579E 

02 

0.57562E 

02 

0.58782E 

02 

0.57403E 

02 

0.58835E 

02 

0.57252E 

02 

0.62612E 

02 

0.57107E 

02 

0.64644E 

02 

0.56985E 

02 

0.64873E 

02 

0.56893E 

02 

0.63293E 

02 

0.56832E 

02 

0.59949E 

02 

C.56797E 

02 

0.57559E 

02 

0.56773E 

02 

0.61782E 

02 

0.56756E 

02 

0.64282E 

02 

0.56757E 

02 

0.64989E 

02 

0.56788E 

02 

0.63884E 

02 

0.56850E 

02 

0.60998E 

02 

0.56940E 

02 

0.56411E 

02 

0.57044E 

02 

0.60831E 

02 

0.57151E 

02 

0.63794E 

81 

0.57271E 

02 

0.64978E 

0.57416E 

02 

0.64351E 

02 

0.57589E 

02 

0.61929E 

0.57781E 

02 

0.57787E 

02 

02 

0.58000E 

02 

97 


HYDE,W.H.   B3   R=0.8  OHMS,  VREF=60.0 


TIME 


VIN 


VOUT 


0.16421E  00 
0.16465E  00 
0.16  5C9E  00 
0.16554E  00 
0.16598E  00 
0.16642E  00 
0.16687E  00 
0.16731E  00 
0.16775E  00 
0.16820E  CO 

0.16864E  00 
0.16908E  00 
0.16953E  00 
0.16997E  00 
0. 17041E  00 
0.17C86E  00 
0.17130E  00 
0.17174E  00 
0.17219E  00 
0.17263E  CO 

0.17307E  00 
0.17352E  00 
0.17396E  00 
0.17441E  00 
0.17485E  00 
C.17529E  00 
0.17574E  00 
0.17618E  00 
0.17662E  00 
0.17707E  00 

0.17751E  00 
0.17795E  00 
0.17840E  00 
0.17884E  00 
0.17928E  00 
0.17973E  00 
0.18017E  00 
0. 18061E  00 
0.18106E  00 
0.18150E  00 

0.18194E  00 
0.18239E  00 
0.18283E  00 
C.18327E  00 
0.18372E  00 
0.18416E  00 
Q.18460E  00 
0.18505E  00 
0.18549E  00 
0.18594E  00 


0.59761E  02 
0.63182E  02 
0.64841E  02 
0.64693E  02 
0.62740E  02 
0.59C37E  02 
0.58576E  02 
0.62448E  02 
0.64578E  02 
0.649C8E  02 

0.63428E  02 
0.60179E  02 
0.55253E  02 
0.48786E  02 
0.40957E  02 
0.31988E  02 
0.22126E  02 
0.11648E  02 
0.84376E  00 
0.0 

0.0 

0.0 

0.0 

O.C 

0.58055E  02 

0.59521E  02 

0.63038E  02 

0.64796E  02 

0.64748E  02 

0.62895E  02 

0.59288E  02 
0.58312E  02 
0.62277E  02 
0.645C7E  02 
0.64937E  02 
0.63557E  02 
0.60405E  02 
0.56988E  02 
0.61396E  02 
0.64091E  02 

0.65000E  02 
0.64096E  02 
0.61405E  02 
0.57002E  02 
0.60394E  02 
0.63551E  02 
0.64936E  02 
0.64510E  02 
0.62285E  02 
0.58324E  02 


0.58213E  02 
0.58431E  02 
0.58665E  02 
0.58922E  02 
0.59200E  02 
0.59489E  02 
0.59775E  02 
0.60055E  02 
0.60343E  02 
0.60645E  02 

0.60962E  02 
0.61286E  02 
C.61604E  02 
0.61894E  02 
0.62130E  02 
0.62282E  02 
0.62313E  02 
0.62185E  02 
0.61878E  02 
0.61535E  02 

0.61194E  02 
0.60855E  02 
0.60517E  02 
C.60182E  02 
0.59848E  02 
0.59517E  02 
0.59187E  02 
0.58863E  02 
0.58562E  02 
0.58286E    02 

0.58030E  02 
0.57780E  02 
0.57534E  02 
0.57308E  02 
0.57111E  02 
0.56945E  02 
0.56806E  02 
0.56681E  02 
0.56562E  02 
0.56462E  02 

0.56392E  02 
0.56356E  02 
0.56350E  02 
0.56363E  02 
0.56383E  02 
0.56418E  02 
0.56481E  02 
0.56575E  02 
0.56700E  02 
0.56846E  02 
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TIME 


VIN 


VOUT 


C.18638E  00 
0.18682E  CO 
0.18727E  00 
0.18771E  00 
0.13815E  00 
0.18860E  00 
0.18904E  00 
0.18948E  00 
0.18993E  00 
0.19037E  00 

0.19081E  00 
0.19126E  00 
0.19170E  00 
0.19214E  00 
0.19259E  00 
0.19303E  00 
0.19347E  00 
0.19392E  00 
0.19436E  00 
0.19480E  00 

0.19525E  00 
0.19569E  00 
0.19613E  00 
0.19658E  00 
0.19702E  00 


0.59276E  02 
0.62888E  02 
0.64746E  02 
0.64799E  02 
D.63045E  02 
0.59533E  02 
0.58042E  02 
0.62101E  02 
0.64430E  02 
0.64961E  02 

0.63681E  02 
0.60626E  02 
0.55880E  02 
0.49575E  02 
0.41889E  02 
0.33035E  02 
0.23259E  02 
0.12835E  02 
0.20527E  01 
0.0 

0.63422E  02 
0.64907E  02 
0.64581E  02 
0.62456E  02 
0.58588E  02 


0.56997E  02 
0.57157E  02 
0. 57338E  02 
0.57546E  02 
0.57781E  02 
0.58035E  02 
0.58293E  02 
0.58550E  02 
0.58820E  02 
C.5  9109E  02 

0. 59419E  02 
0.59744E  02 
0.6007CE  02 
0.60376E  02 
C.60639E  02 
0.60825E  02 
0.60901E  02 
0.60828E  02 
0.60568E  02 
0.60233E  02 

0.59899E  02 
0.59574E  02 
0.59272E  02 
0.58992E  02 
0.58727E    02 
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